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Summary— A review is given of the various possible methods of carrying out batch fractional 


liquid-liquid extraction. Particular attention is paid to the analogies with normal and inverted 
batch fractional distillation and the difference between the two methods of operation is illustrated 
with reference toa specific example. Possible applications of these processes, both in the laboratory 


and in industry, are discussed briefly. 


Résumé— Revue des diverses méthodes de conduire extraction fractionnée liquide-liquide par 


fournées. 


On insiste sur les analogies entre la distillation fractionnée normale ou inversée ; un exemple 
particulier permet dillustrer les différences entre les deux modes opératoires. Bréve discussion 


1. IntTROpUCTION 


Although fractional distillation and liquid-liquid 
extraction are closely related unit operations, 
it is of interest to trace the differences in their 
early development. Thus, although the con- 
tinuous distillation of fermented worts probably 
dates from the development of the Coffey still 
in 1832, most distillation operations were carried 
out batchwise until comparatively recently and 
it is only within the past thirty years or so that 
continuous distillation has become common- 
place. Fractional liquid-liquid extraction, on 
the other hand, has always been considered as a 
continuous or semi-continuous operation, and 
its potentialities as a true batch process do not 
appear to have been recognised. 
The earliest applications of 
extraction were simple single- or multiple- 
contact operations in which charges were 
repeatedly extracted with solvent, similar to those 


liquid-liquid 


des applications possibles du procédé, soit au Laboratoire, soit dans I’ Industrie. 
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carried out by chemists with separating funnels. 
Such methods necessitate a lengthy series of opera- 
tions for the complete separation of mixtures of 
solutes, except in the limit when one of the 
components only is extracted by the solvent. 
In the latter case, no fractionation is required 
and an apparatus similar in principle to the well- 
known Soxhlet extractor is generally used. 


The above methods have given way in more 
recent times to countercurrent extraction, 
analogous to gas absorption, enabling a consider- 
able economy to be effected in solvent usage. 
In 1933, Saat and van Dyck [17] proposed the 
use of extract and raflinate reflux, and in this 
form the process is the complete analogue of 
continuous distillation. When applied to the 
rather rare type of ternary system in which two 
of the three component pairs are partially miscible 
(ef. Fig. 1b), it is possible in principle to obtain 
a complete separation of the components of the 
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feed. Another process employing reflux is that 
in which two or more components are separated 


by fractional distribution between two sub- 


stantially immiscible solvents, although in this 


case the analogy to distillation is less clear since 
produced — internally 
external 


the reflux is generally 
in the extractor, rather than in an 
piece of equipment. This process was used by 
Cornisu et al. [4] for the separation of vitamins, 
and further applications have been described 
by van Dyck [7] and by Scuripes [18]. It is also 
employed in the so-called “ Duo-sol process 
of lubricating oil refining [22]. 

Although, in the past, liquid-liquid extraction 
operations have all been of the continuous type, 
as seen above, a simple laboratory extractor was 
described in 1936 by Cannon and Fenske [3] 
which in operation is strictly analogous to batch 
fractional distillation, although it does not appear 
that this was realised at the time. Recently, 
Weiss and co-workers [15] have deseribed an 
ingenious batch application of the method known 
as “ dissociation extraction ” for the separation 
of mixtures of weak acids or bases. This process 
involves extraction in a column in which a pll 
gradient is maintained by the admission of a strong 
acid at one end and an alkali at the other, thus 
taking advantage of differences in the dissociation 
constants of the substances to be separated. 
The possibility of carrying out such a continuous 
process was first recognised by Janrzen [10], 
who reported results on the separation of mixtures 
of quinoline with isoquinoline and naphthylamine. 
It has been also used by Axe and Bat.ey [2] 
for the separation of complex nitrogenous bases 
present in petroleum, by Scuurze, QUEBEDEAUX 
and Locurr [19] for the separation of n-caproic 
and n-heptylic acids, and by Waker [24] 
for the separation of m- and p- cresols. The 
application of this process to the separation of 
penicillic acids has been described in a patent 
assigned to the Shell Development Co. [20], 
and it has been further discussed by Twice [23] 
and by Ro.ey, Sremer and ZimKin [16]. In 
practice, it is found difficult to maintain sufli- 
ciently accurate control of the pH gradient 
to ensure satisfactory operation under continuous 
process conditions, but the batch method of 


operation of Weiss ef al. shows greater promise 
in this respect. 

As with fractional distillation, it may well 
be the case that there will be a specific need 
from time to time for batch methods of fractional 
liquid-liquid extraction, e.g. when the quantity 
of material to be treated is small, or when a 
plant is required with suflicient flexibility to 
handle charges of widely differing composition. 
Apart from the above twoexamples, batch methods 
of operation do not appear to have received any 
attention in the literature and are very imperfectly 
understood. In this paper, therefore, a description 
is given of the various possible methods of batch 
fractional extraction, pointing out at the same 
time the analogy with known methods of batch 
fractional distillation. It should be noted that 
similar batch methods of operation can in principle 
be applied to ion exchange and adsorption 
processes provided the of adsorbent is 
moved countercurrent to the flow of the fluid 


bed 


stream. 


2. Descrirrion or Processes 
2.1 General 

A large variety of types of equilibria are 
encountered liquid-liquid extraction, and 
systems may conveniently be classified as being 
of type CN, where C refers to the number of 
components present, and N to the number of 
partially-miscible pairs of components. Neglect- 
ing binary (i.e. 2/1) systems, the most important 
of the less complex types of system encountered 
are the 3/1 and 3 2 systems, i.e. ternary systems 
in which one or two pairs of components are 
partially miscible, and 4 1 systems in which 
two solutes are distributed between two sub- 
stantially immiscible solvents. (The latter are 
also occasionally of the 4/2 or 4.3 type). 

Solvent extraction processes, on the other 
hand, are generally classified as single-solvent 
and operations, according as 
one or two solvents are used to effect the separa- 
tion. Thus, in the case of ternary systems a 
single solvent is used to effect the separation of 
the two components of the feed, and one or both 
of the latter may be considered as the second 
solvent, depending whether the system is of the 


double-solvent 
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31 or 32 type. Double-solvent processes, 
on the other hand, are those in which the com- 
ponents of the feed are separated by distribution 
between two substantially immiscible solvents, 
and in the simplest case are of the 4/1 type, 
as stated above. 

The difference between single and double- 
solvent operations may be further clarified by a 
consideration of the phase rule relationships. 
Thus, for a ternary system there are three 
components and three phases (two liquid and one 
vapour), and consequently, according to the 
phase rule, the system possesses two degrees of 
freedom. Hence if the temperature is specified, as is 
usual, the system is defined by fixing a single con- 
centration, which necessarily lies on the miscibility 
boundary. 
there is an additional degree of freedom, and 
consequently two concentrations require to be 
fixed before the system is determined (cf. Fig. 1e). 
The practical result of this is that in double- 
solvent operations the total concentration of the 
solute in the two phases in a continuous extractor 
varies from a maximum at the feed point to 
minima at the column terminals, whereas the 
change in concentration is relatively small with 
single-solvent processes. 


In the case of quaternary systems 


Multicomponent systems become increasingly 
complex in type as the number of components 
is increased, but processes employing these can 
in general be classified as pseudo-single and pseudo- 
double solvent operations. The difference between 
the two may be illustrated by considering the 
case of lubricating oil refining ; thus, the majority 
of processes employ a single solvent such as 
sulphur dioxide, furfural, ete. to extract the 
aromatic from the paraflinic and naphthenic 
components. In the Duo-Sol process, on the 
other hand, the two classes of components are 
separated by fractional distribution between 
propane and a mixture of cresylie acids. 

In the comparison of distillation and extraction 
processes the vapour and liquid phases in the 
former case may be considered to be analogous 
to the extract and raflinate phases in the latter. 
The analogy is completed by the assumption 
that the addition or removal of heat in the 
reboiler or condenser respectively in the case of 


distillation is equivalent to the addition or 
removal of solvent in the case of extraction. 
Thus, the contacting of a charge with a solvent 
in a separating funnel or mixer is analogous 
to the addition of heat to a charge. Both pro- 
cesses are followed by the removal of the new 
phase which is formed, i.e. the extract phase in 
the case of extraction, and the vapour phase in 
the case of distillation. The latter process is 


obviously equilibrium flash vaporisatien. 
Similarly, condensation of vapour, i.e. the 
PLAIT PT 
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Fig. 1. Triangular and equilibrium selectivity diagrams 
for typical systems. 
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removal of heat from saturated vapour to convert 
it into saturated liquid, is analogous to the 
removal of solvent from extract phase, e.g. 
by distillation, to convert it into raflinate phase. 

On the triangular diagram (Figs. la and b), 
the addition of solvent to a mixture of components 
A and B of composition C is equivalent to passing 
along line CS, towards the solvent vertex, as 
far as M, the composition of the resulting two- 
phase mixture. This then separates along a tie- 
line into extract and raflinate phases E and D 
respectively. The removal of solvent from the 
former corresponds to the passage along SE 
on the triangular diagram from E on the extract 
side to F on the raflinate side of the miscibility 
boundary. If further removed, 
analogous to sub-cooling of the condensate, 
ultimately reached. The 


solvent is 


composition Gis 
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extension of the analogy between the two pro- 
cesses to more complex cases will become apparent 
later. 


2.2. Normal batch fractionation 


The process of batch distillation as normally 
carried out using a total condenser is shown 
diagrammatically in Fig. 2(a). In this process, 
the charge in the kettle A is boiled and the 
resulting vapour passed into the fractionating 
column B. The vapour leaving the head of the 
column is condensed in condenser C and _ split, 
part being taken off as product and the remainder 
returned to the head of the column as reflux. 
In principle, the heat removed in the condenser 
can be transferred back to the reboiler by means 
of heat pump D, although in practice this is 
seldom done. 


A MIXER 

(CHARGE 
DISSOLVED 
IN $2) 


EXTRACTION 
COLUMN 


SOLVENT 5, 


SOLVENT 


H MIXER 


(d) DOUBLE SOLVENT EXTN 
CHARGE DISSOLVED IN S; 


(e) DOUBLE SOLVENT EXTN 
CHARGE DISSOLVED IN 


Normal batch fractionational distillation and extraction. 
F indicates position of catch-pot for total reflux stripping. 


2. Broken lines indicate alternative method of double-solvent operation. 
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The analogous process of single-solvent extrac- 
ion is shown in Figs. 2(b) and (c), the former 
for light and the latter for heavy solvents, 
‘onsidering Fig. 2(b), the charge is mixed with 
recirculating solvent in mixer A and passed 
into the base of fractionating column B, the 
indissolved raflinate separating and returning 
to the charge. The extract passes upwards as 
the continuous phase C, is removed at the top 
of the column, and passes to the solvent-removal 
unit, where it is converted into raflinate phase. 
The latter is split, part being taken off as product 
and the remainder refluxed to the top of the 
column, while the solvent is recycled to the charge 
tank. The reflux passes down the column as 
dispersed phase, countercurrent to the extract 
phase, coalesces at the interface and returns 
to the charge tank with the recirculating charge. 
When the solvent density is greater than that 
of the charge the scheme shown in Fig. 1(c) 
is used; this is identical with the foregoing except 
that the column and auxiliaries are inverted. 

Double-solvent batch extraction can be carried 
out in two ways as shown in Figs. 2(d) and (e), 
depending whether the charge is_ dissolved 
initially in solvent S, or S,. Considering the 
former case first and assuming S, to be the 
heavier solvent, the charge dissolved in S, is 
mixed continuously with S, in mixer 4A. The 
mixture then enters the base of the fractionating 
column B where it separates. The S, extract 
passes upwards as continuous phase and on 
leaving the top of the column passes to the S,- 
removal unit. The residue is withdrawn as 
product and solvent S, is recycled to the charge 
tank. Solvent S, enters at the top of the column, 
passes downwards as dispersed phase and coalesces 
at the interface. The S, phase leaving the base of 
the column enters the S,-removal unit E where a 
quantity of S, equivalent to that entering at the 
top of the column is separated and recycled. 
The concentrated S, raflinate then returns to the 
charge tank. When the initial charge is dissolved 
in solvent S, instead of S, the process is identical 
except that the column and auxiliaries are 
inverted as shown in Fig. 2(e). 

The double-solvent processes may be modified 
by splitting the solvent-free S,- or S,-extract, 


taking off a portion as product, dissolving the 
remainder in the other solvent in mixer H and 
refluxing the resulting solution to the column in 
place of the pure solvent feed. This arrangement, 
which is similar to that proposed by AssELIN 
and Comincs [1] in the case of continuous 
extraction, is indicated by broken lines in Figs. 
2(d) and (e). It is likely to possess particular 
advantages in the case of batch extraction since 
it enables the reflux ratio to be altered at will 
without changing the solvent flows. 


2.3. Inverted batch fractionation 


Inverted batch distillation, which is apparently 
a little known process, was first used by HurrmMan 
and Urey [9] for the concentration of the 0' 
isotope by the fractionation of water. It consists 
in distilling the charge in such a way that the 
product is taken off at the base of the column 
and consists initially of the least volatile com- 
ponent of the feed. The remaining components 
are taken off in increasing order of volatility. 
Referring to Fig. 3(a), the charge is fed from 
the charge vessel G into the top of the fractionat- 
ing column B and the liquid leaving the base 
of the column is split, a portion being taken off 
as product. The remainder is vaporized in a 
boiler A of low hold-up capacity, preferably a 
flash boiler, and the vapour is refluxed to the 
base of the column. The overhead vapour is 
condensed in condenser C and the condensate 
is returned to the charge vessel. 

The analogous process for single-solvent systems 
is shown in Figs. 3(b) and (c). Considering first 
the case shown in Fig. 3(b) where the solvent 
is less dense than the charge, the latter is fed 
from the charge tank G into the top of the column 
and passes downwards as dispersed raflinate 
phase, countercurrent to the ascending extract 
phase. The raflinate coalesces at the interface 
and a portion is withdrawn as product. The 
remainder is converted into extract phase by 
dissolving in solvent in mixer A of small hold-up 
capacity and refluxed to the base of the column. 
The extract phase leaving the top of the column 
enters solvent-removal unit C, from which the 
residue returns to the charge and the solvent to 
the mixer. The arrangement shown in Fig. 3(c) 
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Fig. 3. Inverted batch fractionation. 


Notes: 1. F indicates position of catch-pot for total reflux stripping. 


2. Broken lines indicate alternative method of double-solvent operation. 


is used when the solvent is denser than the feed, 
and is identical except that the equipment is 
inverted. 

The arrangement used for double-solvent 
processes is shown diagrammatically in Figs. 3(d) 
and (e), depending whether the charge is dissolved 
initially in solvents S, or S,. Considering Fig. 3(d) 
first, the charge in vessel G, dissolved in solvent 
S,, is fed to the top of column B and passes down 
as dispersed phase countercurrent to an upflow 
of solvent S,. The S, phase leaving the base of 
the column enters the S,-removal unit E and the 
residue is removed as product. Solvent S, 
enters the base of the column and leaves at the 
top, entering the S,-removal unit C. The S, 


solvent which is removed is recycled to the base 
of the column, while the residue is mixed with 
recycled solvent S, and returned to the charge. 
In the alternative scheme indicated by broken 
lines, the product stream, after removal of S,, 
is split and a portion is dissolved in S, and refluxed 
to the column. The arrangement shown in Fig. 
3(e) is identical except that the charge is initially 
dissolved in solvent S,. 


2.4 Total reflux stripping 


When it is desired to separate a very small 
amount of more-volatile component from a 
charge a method of batch distillation known as 
“total reflux stripping” is sometimes applied, 
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COLUMN 


For this purpose a small catchpot is inserted in 
the reflux line at F, Fig. 2(a). The distillation 
is carried out under total reflux conditions, without 
withdrawal of product, until equilibrium is 
reached, so that the catchpot contains the 
highest possible concentration of low-boiling 
component. The contents of the catchpot are 
then withdrawn as product, and the process is 
repeated. 

It is apparent that total reflux stripping can 
be applied to single-solvent extraction in an 
indentical manner to that used for batch distilla- 
tion, by inserting a catchpot at F in Figs. 2(b) 
and (c). It can also be used in conjunction with 
double-solvent extraction, provided the alterna- 
tive method of operation indicated by broken 
lines in Figs. 2(d) and (e) is used, and the catchpot 
is inserted at F as before. 

Total reflux stripping can also be used in 
conjunction with the inverted batch fractionation 
methods shown in Fig. 3. Thus, in the distillation 
case the catchpot is inserted below the column at 
F in the liquid withdrawal line. The position 
of the catchpot in the case of single- and double- 
solvent extraction is similarly indicated at F 
in Figs. 3 (d) and (e). 


2.5 Dissociation extraction 


The process of dissociation extraction as applied 
to the separation of a mixture of weak acids 
is shown in Fig. 4. In the continuous form 
of this process shown in Fig. 4(a) the aqueous 
solution of acids is fed into the mid-point of the 
extraction column and aqueous caustic soda 
is fed at the top. Similarly, organic solvent 
together with a mineral acid is fed at the bottom 
of the column and the quantities of acid and 
and alkali are adjusted so that a pH gradient is 
maintained throughout the length of the extractor. 
As a result, the acid which has the lower dissocia- 
tion constant is extracted preferentially by the 
organic solvent, while that with the higher 
dissociation constant passes into the aqueous 
product at the bottom of the column. 

The batch form of this process, as described 
by Wess et al. [15], is shown in Fig. 4(b); 
the process in this form corresponds to that shown 
in Fig. 2 (ec). The charge, dissolved in organic 
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Fig. 4. Dissociation extraction. 


solvent, (hereafter referred to as “solvent ”’) 
is contained in the charge vessel and is contacted 
with aqueous caustic soda. The mixture passes into 
the top of the fractionating column and the organic 
layer coalesces at the interface. The aqueous 
layer passes down the column as continuous 
phase and on leaving enters the water-removal 
unit, corresponding to the S,-removal unit in 
Fig. 2(e), which in the present case is an extraction 
column with a feed of solvent and acid entering 
at the base. The stripped water together with 
sodium sulphate leaves at the bottom of this 
column, while a solution of the acids in solvent 
leaves at the top. This stream is split, part 
being taken off as product and the remainder 
returned to the fractionating column. The 
solvent phase passes upwards as the dispersed 
phase and coalesces at the interface where it 
mixes with the circulating solvent solution of the 
charge. The solvent layer is taken off at the top 


WATER +Na»SO,4 H7S04 FEED 


of the column and enters a_ solvent-removal 
unit which, in this case, is a continuous still. 
The residue, consisting of a concentrated solution 
of the acid in solvent, is returned to the charge 
tank, while the solvent recycles to the base of 
the water-removal column. This process is 
basically similar to that shown in Fig. 2(e), 
but it is obvious that a similar process could be 
devised in which the charge is disolved initially 
in the aqueous solution ; inverted batch methods 
similar to those shown in Figs. 3(d) and (e) 
could also be used. 


2.6 Acilliary equipment required 


In batch distillation three main items of equip- 
ment are required, namely the charge tank and 
reboiler, which are generally but not necessarily 
combined, the fractionating column, and the 
condenser. Similarly, in solvent extraction three 
items are required for single-solvent processes, 
the charge tank and solvent mixer, which again 
may or may not be combined, the fractionating 
column, and the solvent-removal unit; in the 
ease of double-solvent processes a second solvent- 
removal unit is required. 

The fractionating column used in batch solvent 
extraction operations may be of any standard 
type, e.g. spray, packed, mechanical or pulsed ; 
mixer-settlers may also be used, provided the hold- 
up is not excessive. A number of methods of 
extracting the charge with solvent in normal 
batch operation is shown in Fig. 5. The scheme 
shown in Fig. 5(a) is the simplest and consists 
in using the base of the column as the charge 
vessel, a spray of solvent being introduced at the 
base. This system was used in the laboratory 
extractor described by Cannon and Fenske [3]. 
Similar methods in which a mechanical agitator 
is used to improve the contacting are shown in 
Figs. 5(b) and (c); the latter has the advantage 
that the initial charge can be reduced down to a 
very small volume. More elaborate systems in 
which the charge is contained in a_ separate 
vessel are shown in Figs. 5(d) to (f); in the first 
two, the solvent is mixed with the charge in the 
charge vessel and the mixture is introduced 
to the separating section of the column, either by 
gravity or by means of a pump. All the foregoing 
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(¢) (e) (f) 
Fig. 5. Methods of introducing solvent (ef. Figs. 2b. and 
2d). 


methods are analogous in the case of batch 
distillation to the use of a kettle with an internal 
steam coil. The system shown in Fig. 4(f), 
on the other hand, corresponds to the use of a 
separate charge tank and external calandria and 
consists in circulating the charge and solvent 
simultaneously by means of a centrifugal mixing 
pump into the base of the column, where phase 
separation takes place. 

In the case of inverted batch distillation it is 
necessary to maintain the hold-up in the reboiler 
as low as possible, and for this reason a flash 
or film evaporator is normally used. Similarly, 
with inverted batch extraction it is necessary 
to reduce the volume in the mixer to a minimum 
and for this reason a centrifugal pump would 
probably prove most suitable. 

The solvent removal unit can take a number 
of forms, depending upon the type of system 
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used and the physical properties of the components. 
In general, continuous distillation columns, 
liquid-liquid extraction units or simple evapora- 
tors are likely to be the most useful ; with double- 
solvent operations the two solvent removal units 
may sometimes be of different types. 

With  single-solvent systems the solvent 
removal unit would generally consist of a con- 
tinuous distillation column in which the solvent 
is taken off either as bottom or top product, 
depending upon its volatility relative to that of 
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Fig. 6. Types of solvent removal unit (cf. also Fig. 4b 
for use of simple extraction column). 
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the solute (cf. Fig. 6a). A similar method can be 
used for double-solvent operations with volatile 
solutes, but with non-volatile solutes a simple 
evaporator would be used, as shown in Fig. 6(b). 
When it is necessary to transfer the solute from 
one solvent to the other, as for example when the 
alternative method of providing reflux indicated 
by broken lines in Figs. 2 and 3 is used, or when 
returning the extract to the charge tank in 
inverted double-solvent operation (Figs. 3d and 
3e), an extraction column can be used, as shown 
in Fig. 4(b) for the case of dissociation extraction. 
Alternatively, with non-volatile solutes a 
combined evaporator and extraction column 
can be used, as suggested by AsseLIN and 
ComminGs [1]; an example of the latter arrange- 
ment in which a climbing film evaporator is used 
to obviate the need for a pump is shown in 
Fig. 6(c). This type of evaporator has the addi- 
tional advantage that the hold-up is very low, 
and this is an essential requirement of all types 
of solvent removal unit since a large hold-up 
will reduce the sharpness of separation of the 
components, as with batch distillation. 


3. Discussion 


Considering first single-solvent operations, there 
is generally little point, as mentioned earlier, 
in using reflux with systems of type 3/1 since the 
degree of enrichment of the extract product 
which can be obtained is limited by the plait point 
composition. Consequently the use of the pro- 
cesses shown in Figs. 2(b) and (c) and 3(b) and (ec) 
is likely to be restricted in practice to the rather 
rare type 3 2 system. The system of this type 
which has most generally been considered is 
aniline-cyclohexane — n-heptane, but it is doubtful 
whether this has ever been employed in large 
scale practice. Amongst other type 3/2 systems, 
few, of those investigated up to the present are 
likely to have any industrial application, with 
the notable exception of the liquid propane- 
fatty acid systems investigated by H1rxson 
and co-workers [6, 8]. 

Double-solvent processes, on the other hand, 
are potentially of considerably greater interest, 
and their use has been described for the separation 
of vitamins [4], of azeotrope-forming mixtures 


such as ethanol and methyl ethyl ketone [18], 
and of mixtures of isomers, e.g. o- and p- 
chloronitrobenzene [18]. They are also becoming 
of rapidly increasing importance for inorganic 
separations, e.g. in the separation of rare earths 
by extraction with tributyl phosphate solutions 
[25]. Although these processes were originally 
described in continuous form, they can easily be 
translated into the batch forms shown at (d) 
and (e) in Figs. 2 and 3. Processes of this type 
ean be carried out in the laboratory on a very 
small scale, particularly when using small! diameter 
mechanical columns of high efficiency, and they 
may well prove to be of comparable value to the 
Craig machine [5] for small-scale analytical 
separations. On a somewhat larger scale they 
could be used with advantage for investigating 
the possibility of carrying out full-scale separa- 
tions by solvent extraction, and in this respect 
they bear the same relationship to continuous 
extraction as laboratory batch distillation does 
to full scale continuous distillation. On a still 
larger scale they should prove of considerable 
value for the small-scale production of valuable 
materials in, e.g., the pharmaceutical industry, 
and for this purpose the flexibility of batch 
operation is particularly advantageous since a 
standard piece of equipment could be used for 
a range of products. 

The application of the four processes shown 
in Figs. 2(d) and (e), and 3(d) and (e), can best be 
illustrated with respect to a specific separation 
e.g. the separation of o- and p-chloronitrobenzene 
by fractional distribution between a heptane 
mixture and aqueous methanol [18]. Two of 
these processes require the charge to be dissolved 
initially in the methanol phase, and two in the 
hydrocarbon phase. Considering the former 
case first, the process shown in Fig. 2(d) consists 
in dissolving the chloronitrobenzene mixture in 
85°, methanol, contacting with heptane in 
mixer A and feeding to the base of the column. 
The heptane phase leaving the top of the column, 
which is enriched in the para isomer, passes to 
the solvent-removal unit, and part or all of 
the residue is taken off as product, depending 
upon the method used to provide the reflux. 
As the fractionation continues the concentration 
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of the ortho isomer in the product will increase, 
and the reflux ratio is consequently increased 
gradually to maintain a high degree of separation; 
however, a point is reached where this is no 
longer practicable and the product is “ cut” 
to the intermediate fraction receiver. The 
process is then continued until substantially the 
whole of the more soluble para isomer is extracted 
and the residue, consisting of a solution of ortho 
isomer in aqueous methanol, is discharged as less- 
soluble product. Alternatively, the extraction may 
be continued until the whole of the ortho isomer is 
separated, leaving a residue of aqueous methanol. 

The same charge may also be treated by 
inverted fractionation using the scheme shown 
in Fig. 3(d). In this process the product consists 
initially of a solution of the ortho isomer in 
aqueous methanol, leaving a residue of the 
para dissolved in the same solvent. Thus, 
the charge is fed to the top of the column and is 
extracted by an ascending stream of heptane. 
The methanol phase leaving the base of the 
column passes to a solvent removal unit, and part 
or all of the residue is taken off as product. 
Fresh heptane, in which the remainder of the 
residue is dissolved if the reflux is provided 
externally, is fed to the column base. The 
heptane phase leaving the top of the column 
enters the second solvent-removal unit of the 
type shown in Fig. 6(c), and the residue passes 
as a methanol solution back to the charge tank. 
As the extraction progresses the concentration 
of para isomer in the product tends to rise, and 
the reflux ratio is gradually increased to prevent 
this. Ultimately the para content becomes excess- 
ive and the product is “* cut ” to the intermediate 
fraction receiver and the extraction continued 
until it is free from ortho. The remainder of the 
ortho isomer is then recovered, either directly 
from the methanol solution in the charge tank, or 
by continuing the extraction as before. 

For the remaining two processes the charge 
would be dissolved initially in the heptane phase 
This would then be processed as before by normal 
or inverted fractionation, when the initial product 
would consist either of pure ortho or pure para 
isomer respectively. 

The decision as to which of the above four 
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methods of operation should be used in any 
particular case depends upon the nature of the 
problem. Where possible, the desired component 
of the mixture should be taken off as initial 
product, particularly if it is present in smaller 
amount, since the separation would normally 
be sharper under these circumstances. Having 
decided which component is to be taken off 
initially there is then a choice between using 
normal fractionation with the initial charge 
dissolved in one of the solvents or inverted 
fractionation with the charge dissolved in the 
other solvent. In some cases it may be found 
that inverted operation is the more convenient from 
the practical point of view, since then the mixer is 
not required. In other cases, the type of operation 
may be dictated by the previous history of the 
charge. This is particularly true of inorganic 
solutes where it is normally preferable to start 
with an aqueous solution. It is apparent, there- 
fore, that at the present time definite rules 
cannot be given as to which process to adopt, 
but the position will no doubt be clarified as 
experience is gained with operating techniques. 

With batch processes in general, intermediate 
fractions consisting of mixtures of two components 
are always obtained between the main component 
fractions and these would normally be kept to a 
minimum by increasing the reflux ratio in steps 


as the purity of the main fraction decreases as 
described above. When processing repeated 
batches of similar composition, the intermediate 
fractions would be recycled to the next charge ; 
by this means the fresh charge is split completely 
into its components and intermediate fractions 
of the same size and composition as those recycled 
to the charge vessel are obtained. If the holdup 
in the column is negligibly small in relation to the 
size of the charge, batch fractionation curves may 
be calculated by graphical integration of the 
Rayleigh equation, as described by Smoker and 
Rose [21] for the case of batch distillation, 
making use of the curve relating the product to 
the charge composition. The latter may be 
obtained for single-solvent systems by the pro- 
cedure described by Matoney and Scuvusert [14], 
and for double-solvent systems by the use of 
methods described by Scuripe. [18] and by 
KLINKENBERG and co-workers [11-13]. 


4. ACKNOWLEDGMENTS 


Especial thanks are due to Messrs. Wetss 
and Swinton of the C.S.1.R.O., Melbourne, 
for demonstrating their method of batch dissocia- 
tion extraction and for stimulating the writer 
to continue his study of batch processes in general. 
Thanks are also due to Mr. A. S. Wurre for his 
interest and encouragement in this work. 


REFERENCES 
{1] Assen, G. F. and Comines, W.; Ind. Eng. Chem 1950 42 1198. 
(2] Axe, W. N. and Barwey, J. R.; J. Amer. Chem. Soc. 1939 61 2609; ef. Axe, W. N., ibid. 1017. 
[3] Cannon, M. R. and Fenske, M. R. ; Ind. Eng. Chem. 1936 28 1035. 
[4] Connisu, R. E., Ancuracp, R. C., Murpny, E. A. and Evans, H. M.; Ind Eng. Chem. 1934 26 397. 
(5) Crase, L. C. ; J. Biol. Chem, 1943 150 33 ; ibid. 1944 155 519; Crate et al. ; ibid. 1945 161 321; 1947 168 665. 
(6) Drew, D. A. and Hixson, A. N.; Trans. Amer. Inst. Chem. Eng. 1945 40 65. 


[7] Dyex, W. J. D. van. ; Pet. Zeit. 1982 28 (8) 7. 


[8] Hixson, A. W. and Bocke mann, J. B. ; Trans. Amer. Inst. Chem. Eng. 1942 38 891. 
[9] Hurrman, J. R. and Urey, H.C. ; Ind. Eng. Chem. 1937 29 531. 
(10] Janrzen, E.; “ Das fractionerte Destillieren und das fractionerte Verteilen,” Dechema Monographie, Band 5, 


No. 48, pp. 81-117. (Berlin, Verlag Chemie, 1932). 
Kuwkenperc, A. ; Chem. Eng. Sci. 1951 1 86. 


(12) KuwKensere, A., Lauwerter, H. A., and Reman, G. H.; Chem. Eng. Sei. 1951 1 93. 


[13] Kumwkenperc, A.; Ind. Eng. Chem. 1958 45 653. 


[14] Maroney, J. O. and Scuuperrt, A. E. ; Trans. Amer. Chem. Eng. 1940 36 741. 
{15] Minty, D. W., McNer, R., Ross, M., Swrnron, E. A. and Weiss, D. E., Australian J. App. Sei. 1953 4 530. 


199 


H.R. C. Prarr: Batch fractional liquid-liquid extraction 


[16] Roney, D., Sremer, H. and Zim«ry, E.; J. Soe. Chem. Ind. 1946 65 237. 

[17] Saat, R. N. J. and vaw Dyck, W. J. D.; Proce. World Pet. Cong. (Lond.) 1933 11 352. 

[18] Scnermer, KE. G.; Chem. Eng. Prog. 1948 44 681; ibid. 771. 

{19] Scnurze, H. G., Quenepraux, W. A. and Locure, H. L. ; Ind. Eng. Chem. (Anal.) 1988 10 675. 

[20] Shell Development Co. ; B.P. 599, 854 (1948). 

[21] Smoker, E. H. and Rose A.; Trans. Amer. Inst. Chem. Eng. 1940 36 285. 

[22] Torrie, M., and Minver, M. B.; Proc. Amer. Pet. Inst. 1933 14 (II), 85. 

[23] Twies, G. H.; Nature 1949 163 1006, 

Waker, C. A.; Ind. Eng. Chem. 1950 42 1226, 

[25] Weaver, B., Kaprecmann, F. A., and Torr, A. C. : ORNL-1408 (22/10/52) ; Torr, A. C. ; ORNL-1409 (26 /9 /52). 


VOL. 
3 


1954 


Chemical Engineering Science, 1954, Vol. 3, pp. 201 to 208. Pergamon Press Ltd. 


Combustion rate of coal shale 
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Instituut voor Industriéle Scheikunde, Universiteit te Leuven 


(Recewed 9 July, 1954) 


Summary-—The combustivity of coal shale and other materials containing combustible carbon 
disseminated in the non-combustible matter, is governed by the possibility of access of oxygen 
to the burning surface through the porous layer of mineral matter left by the combustion. In 
order to express the combustion rate, a relation is established, one of the terms being the 
diffusivity through the residual porous material ; experimental values were obtained for the latter. 


Résumé— La combustibilité des schistes de lavoirs et d'autres matiéres contenant du carbone 
combustible disséminé dans une substance incombustible, est déterminée par la possibilité 
d'accés de loxygéne a la surface en combustion a travers la couche poreuse de matiéres minérales 
laissée par la combustion. Pour exprimer la progression de cette combustion, une relation 


est établie dont un des termes est la diffusivité & travers la matiére poreuse résiduelle ; des valeurs 


expérimentales ont été obtenues pour cette diffusivisé. 


The increasing necessity of improving the overall 
economic pattern of the coal industry has led 
to several suggestions, one of them consisting 
in the combustion of coal shale, the waste product 
of coal cleaning. 

Coal shale contains some 10 to 20°, imbedded 
combustible material. The combustion of this 
poor fuel is a particular problem ; it is obvious 
that the possibility of combustion is governed 
by the accessibility of the imbedded combustible 
material to oxygen : in other words, the porosity 
of the material is a governing factor. As a 
direct consequence, the combustion temperature 
must be kept low enough to avoid sintering. 

A very similar problem has been encountered 
in the shale oil industry, where a residue is left 
after pyrolysis of the oil shale, commonly 
described as “ shale coke,’ and containing some 
5% combustible carbon distributed throughout 
the stony material. In a _ well-known recent 
process, boiler tubes are inserted in the burning 
bed of shale coke, in order to keep the main 
temperature below the sintering temperature 
[3] ; the applicability of a similar process to coal 
shale combustion is under investigation [4]. 
Design and operation of such kilns would be 
facilitated by the knowledge of combustion time 
at various conditions ; the purpose of the present 


study is to allow the prediction of the progress 
of combustion for different kinds of coal shales. 

The rate of reactions accompanied by the 
formation of a porous layer has been investigated 
so far with special attention to the reduction of 
ores [1] [11]. More recently, a general treatment 
of reactions involving diffusion and simultaneous 
gas flow has appeared [5]. However, the com- 
bustion rate of coal shale and similar materials 
presents some particularities seemingly justifying 
a separate treatment. 

The rate of combustion of combustible material 
ig controlled by two factors: the rate of the 
chemical reaction at the combustion surface, 
and the rate of diffusion of oxygen and combustion 
gases respectively to and from the reacting surface. 

The combustion of carbon, expressed in terms 
of oxygen consumption at the carbon surface, 
is a first order process at temperatures below 
1500°C, and can be expressed by : 


N, =k. (1) 
where: Ny = molar rate of consumption of 


oxygen per unit area of surface. 


T, = temperature of carbon surface. 


E = energy of activation, depending 
on the sort of carbon 
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Po, = partial pressure oxygen at carbon 
surface 


k =a constant 


The second factor is gas diffusion. In the case 
of a material burning as a whole, the rate of 
access of oxygen to the carbon surface is defined 
by the diffusion rate through a stagnant gas 
film contiguous to the carbon surface. However, 
in the present case, the main resistance offered 
to the passage of gases is to be expected as coming 
from the porous layer of non-combustible material 
surrounding the burning carbon, It seems reason- 
able to admit that the diffusional resistance of 
the gas film is negligible in comparison to the 
diffusional resistance of the porous layer: as a 
matter of fact, the validity of this assumption 
will be confirmed by the measurements described 
in the present paper. 

Following the classical treatment for gas diffusion, 
but bearing in mind the particular circumstances 
of the case here considered, it is easily shown that 
the molar rate of oxygen diffusion can be expressed 
as : 

RT py 


L representing no longer the thickness of the 
diffusion film, but the variable thickness of the 
layer of inert material, D being the diffusivity 
under present conditions; P is total pressure, 
Pog 2nd po, partial pressures of oxygen respectively 
in gas phase and at reaction surface, p, the usual 
“ film pressure ” becoming a “ layer pressure.” 

Considering equimolar counter diffusion of 
O, and CO,, neglecting other effects such as 
moisture and volatile matters diffusing outwards, 
the foregoing expression reduces to : 


N 


» D 
N, = RT (Po — Po) (2) 
Combining (1) and (2), eliminating pp, : 


N Pog (3) 


This equation has the well-known form 


driving force 
~ diffusional resistance + chemical resistance 
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Several authors, among others Horren: and 
co-workers [7], recognized that at temperatures 
around 1,100°K and _ higher, the chemical 
resistance is negligible in comparison to the 
diffusional resistance; this applies even more 
firmly to the present case, since the diffusional 
resistance of the layer of inert material is obviously 
much higher than the diffusional resistance of the 
gas film. Then (3) becomes : 


D p 

writing p for po, since there is no longer 
ambiguity. Introducing the specific combustion 
rate C, defined as the weight of carbon burning 
per unit time and surface, it is seen that 

Cc = 12 Ny 


assuming that all the carbon is burnt to CO,. 
The validity of this assumption may be questioned : 
it has long been recognized that at least part 
of the carbon burns primarily to CO, and more 
recently it has been suggested from different 
sides [1, 5, 6, 8) that the primary combustion 
of carbon leads substantially to CO, But it must 
be taken into account that the combustion of the 
primary product CO is a chain reaction and 
proceeds at a sufliciently high rate to mask almost 
entirely the two-step combustion mechanism, 
if no precautions are taken, Consequently, 
it seems permissible to consider only the overall 
reaction, which amounts to admitting that the 
combustion of the CO, following immediately 
the formation of CO, is achieved in the immediate 
vicinity of the carbon surface, so that only CO, 
is involved in the diffusion process. We may 
write accordingly, 


12D p 

C, = (4) 
where L, the thickness of the layer of inert 
material, varies with combustion time, and must 
be expressed in terms of this variable. 

Consider therefore a lump of shale, of perfect 
spherical shape, with radius r,; suppose that 
after time @ the combustion zone is delimited 
by a sphere with radius r, containing the remainder 
of the carbon; the material is assumed to be 
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homogeneous. Equation (4) is not readily 
applicable, since it has been established on the 
assumption of an infinite plane surface. Allowance 
for the fact that the surfaces under consideration 
are spherical is easily introduced, reasoning as 
Horre. did in a similar case. It amounts to 
expressing both sides of (4) in terms of the relevant 
surfaces. For the left-hand side, this surface is 
the area of the sphere enclosing the unburned 
carbon; for the right-hand side, the relevant 
surface is the geometric mean of the outer surface 
of the lump and the surface of the sphere 
containing the unburned carbon. Thus : 
1 
C,. 427? = RT 'L 
and, observing that L = r, —r and simplifying 


[4x (r + 40 


5 
RT —ryr (5) 


The rate of combustion of a sphere with density 
p and containing c,°,, of carbon can be expressed 
by its rate of change of volume times its density 
times its carbon content : 


dr/dd 


but also by its surface times its specific com- 
bustivity : 
. C,. 


Equating the two expressions and simplifying : 


Introducing the values of C, from (5) and 
rearranging : 
(% —r)rdr 1200 DP ag. 
To RT pc 
On integration this yields : 
(ry — 1)? + 2r) 1200 D p 
6 RT peg 


(6) 


This relation between r and @ is not readily 
applicable, since the experimental information 
on the rate of combustion will usually be expressed 
in terms of the residual carbon content rather 


than by values of r. In order to set up the 
necessary relationship, let ¢ be the carbon content 
after time @. The quantity of carbon left in the 
shale can be expressed either by the total weight 
at time @ times c/100: 


ec 4 
100 3770" Pe 

p, being the density at time @, or by the weight 

of the unburned sphere times its carbon content 

100 


Equating and simplifying : 


Call c’ the carbon content of the original material; 
it is seen that ec’ can be deduced from ¢ by the 
relation 
— [100 — + 

100 — (e + C”) 


C” being the CO, content due to carbonates. 
Then by definition 


..p=C. 
and 
or r/t, = (c’/e,)3, 


rearranging this relation between r and the 
experimental values ec’ and ¢, : 


Introducing the expressions for r/r, and (r, — r) 
in (6) gives : 

[1—(e’ /eg)4]? . [1 +2 = RE pear (7) 
The kinetic theory indicates that diffusivity 
varies with temperature approximately according 
to the relation : 


({T 
D,  \273 
where n may have values between 1-5 and 2-0. 


GILLILAND [2] recommends the use of an equation 
in which n is 1-5. Assuming that the relation 
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holds in the present case, and introducing 
D, in place of D, 


[1 — . [1 + 2 = 
(8) 


Equation (8) is the required relation, expressing 
the carbon content c’ after time @ as a function 
of the time of combustion, the partial pressure 
of oxygen in the gas phase, the combustion 
temperature, and constants depending on the 
material considered. Conversely, it can be used 
to calculate the time necessary to reach a given 
combustion ratio; more particularly, the time 
necessary for complete combustion is given by 


. 
19-4 p. 7°. De 


EXPERIMENTAL 


Experimental determinations of combustion rates 
were carried out for five samples of coal shale 
“ coke,” some of the work 
For each 


and one oil shale 
being described in another paper [4]. 
determination, a sample of about 
burned in an electrically heated furnace, in an 
air stream. After combustion, the remaining 
carbon content was determined by combustion 
with oxygen and measurement of the CO,, 
the temperature being carefully maintained 
below the sintering temperature; before the 
analytical combustion, but after weighing, the 
carbonates were broken down by acid. 

Table I gives the properties of the materials 
tested, whilst Table I] summarizes the results. 


75g was 


Table I. 
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900° 
COAL SHALE A © e 
COAL SHALE D 4 a 
COAL SHALE © 


CARBON CONTENT,c —e 


Curves were drawn relating carbon content to 
combustion time at given temperatures, and from 
these curves smoothed values were deduced for 
further calculations, the first step being the 
calculation of the carbon content reduced to the 
original weight ¢c’ by the relation 

100 — + 


100 —(e + 


where ¢, and ¢ are combustible carbon contents, 
CC,’ and C” the CO, in carbonates, respectively 


Properties of shales. 


Material 


Coal shale A 
” B 

D 


\ 
\ 
\ 
| 
\ 
| 
| 
\ \ 
\ \ 
\ \ 
AS — --$———— 
\ \ | 
\ 
a \ | 
fa = ~ 
~ 
~ ~ 3 
= 
1954 
60 120 min 180 
TIME 
Fig. 1. 
FS Density Av. equiv. Carbon content COg in 
p radius total combust. carbonates 
gicn® To cm Co % fo % Cs % 
2-47 0-59 9-79 9-03 2-79 
2-25 0-635 13-08 12-39 2-53 
2-49 0-675 13-68 13-12 2-05 
2-40 0-46 14°80 13-73 3-92 
» 8 2-17 0-49 | 22-09 21-58 1-87 
Oil shale “ coke ” | 2-03 0-61 | 702 4°25 | 10-15 
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Plate I. Partially burned oil shale * coke.” 


(a) Section parallel to stratification. 


(b) Section perpendicular to stratification. 


» 
AK 
4 
(a) 
(b) 
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Table 11. Combustion rate. 


Temperature 


before and after combustion. Uusally, C”’ 
can be neglected, except for oil shale “ coke ” 
burned at 700°C ; the values listed in Table IV 
were used for the latter. 


Discussion 
The value to be used for the temperature term 
in equation (8) is some mean value between the 
temperature of the outer shell and the temperature 


of the carbon surface. It is believed that the 
difference between those two temperatures is 
not very large, and that by taking the temperature 
of the outer shell, which is readily accessible for 
measurement (optically for instance), no large 
error will be introduced, particularly since the 
0-5 exponent minimizes the effect of this error : 
an error of 100° in 1,000°K will lead to a difference 
of about 3°, in the combustion time, which is well 
below the magnitude of accuracy for other 
experimental values. 

The pressure term p is the partial pressure of 
oxygen in the gaseous phase surrounding the 
burning shale. With infinite excess of air, 
p is constant at 0-21 P, and operating at atmos- 
pheric pressure, p = 0-21 atm. With a limited 
excess of air, defined as (1 + ¢) times the theore- 
tical amount, it is easily seen that 

p = 021 + 
l+e 
when the carbon content has changed from 
to ec’. 

The present treatment is based on the assump- 
tion of a spherical particle; clearly the actual 
shape of the shale lumps is quite different. This 
difficulty can be removed by considering r, 
as the equivalent radius defined by 
on lie’ 
P, being the mean weight of a lump of shale. 

In addition, a constant value for r, throughout 
the combustion process has been implicitly 
assumed ; in fact some shales do split during 
combustion, altering the equivalent radius. An 
uncertainty is introduced in this way, since it 
cannot be foreseen if a given shale will split or 
not ; certainly the amount of splitting cannot be 
predicted. In our experimental method, the 
shales were brought abruptly to high tempera- 
tures, in contrast with ordinary practice, so that 
the effect of splitting will be less pronounced in the 
latter case as compared with our experiments. 

It has also been assumed that combustion 
proceeds uniformly from the outside to the 
centre of an individual lump. Somewhat un- 


expectedly, this assumption is in good agreement 
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‘ime content c 
(min.) % 
Coal shale A 200 60 2-93 
120 1-50 
180 0-56 
900 60 2-03 
” 120 0-58 
” 180 0-51 
Coal shale B 800 30 6-87 
” 60 5-05 
” 120 2-60 
900 30 6-61 
60 4-97 
” 120 2-20 
"eo Coal shale C 800 60 3-81 
” 120 0-71 
4 900 60 3-65 
” 120 0-62 
Coal shale D 800 60 5-15 
” 120 2-81 
” 180 0-83 
900 60 4-74 
" 120 1-74 
” 180 0-88 
Coal shale E 800 60 11-21 
” 120 4:10 
180 0-84 
900 60 7-04 
” 120 2-18 
o 180 0-90 
Oilshale“coke” 7000 80 
” 60 2-03 
” | 120 1-32 
800 | 30 | 0-90 
60 0-24 
|_| 
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Table 3. Experimental values for diffusivity. 


Combustible 
Material Combustion Combustion carbon Diffusivity 
temperature time a content b Dy 
T K O sec. (reduced) em sec. 
c % 
Coal shale A 1073 1s = 10 30-04 48 0-095 3-1 x 10% 
36 =< 10° 61-88 2-7 0-255 41 x 10° 
123-75 12 O-485 39 10°3 
1173 1s = 108 32-32 41 0-136 42 x 10°3 
3-6 = 108 64-64 20 0-343 5-3 x 10°43 
7-2 x 10° 129-28 0-640 49 104 
Coal shale B 1073 1s = 10° 21-40 71 0-076 36 « 10-3 
3-6 108 42-80 45 0-200 46 108 
72 x 108 85-50 2 0-402 47 x 10°3 
1173 1s 108 22-38 6-8 0-087 59 = 10°8 
3-6 10% 44-70 0-205 46 = 108 
72x 108 89-50 0-446 5-0 = 10°73 VOL. 
Coal shale C 107% 10 16-15 O-122 75 1954 
36 10° 32-30 B35 0-293 91 10°43 
7-2 108 64-60 0-709 10-9 
1173 18 108 16-88 6-0 O-135 80 
36 108 33-75 3-3 0-310 92 x 104 
72 =x 108 67-50 O4 O-763 11-3 x 10°74 
Coal shale D 1073 18 = 108 34-30 71 0-102 3-0 x 10°43 
3-6 x 108 68-60 45 0-228 32 x 107-3 
72x 10 137-20 21 0-448 33 x 10-4 
1173 8 37-40 65 O-125 23 x 10°43 
3-6 108 T4-80 41 0-257 3-4 x 10°43 
72 x 108 149-60 18 0-488 3-3 x 1073 
Coal shale E 1073 18 « 108 20-96 14:3 0-045 21 x 10-4 
36 x 108 41-92 o1 0-155 37 x 10-4 
72 108 3-2 0-455 54 10-4 
1173 = 108 1S 0-004 4:3 x 10-4 
36 « 10% 43-68 0-289 66 10% 
72 x 103 87-36 18 0-590 67 x 
Oil shale coke 97: 18 = 10° 715 2-54 0-069 0-96 = 10% 
3-6 10° 143 1-92 0-137 0-96 
108 286 1-23 0-265 0-93 10-4 
1073 10 73 O79 0-395 5-42 x 10°4 
3-6 = 108 145 0-21 0-696 +380 10-4 
b = [1 — . [1 + 2 


with the actual facts; the progress of the com- perpendicular to stratification. This is shown by 
bustion zone is nearly the same in the direction the illustrations. 
of the stratification planes as in directions There remains the term of diffusivity through 
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the porous layer. From the experiments, values 
were obtained for D,, as summarized in Table III. 
Considering the assumptions introduced, the 
values obtained for D are remarkably consistent 
for coal shale D, reasonably good for coal shales 
A and B, especially if the values after } hour com- 
bustion are discarded, considering the lower 
accuracy of their determination. For shales ( 
and E, a more systematic increase of D, with 
combustion time and temperature seems to occur. 
This can be reconciled to the observation that C 
and E showed a good deal of splitting, whilst 
this effect was less pronounced for A and B, 
and practically absent for D. Thus the variation 
in equivalent diameter during combustion as a 
result of splitting may be responsible for the 
observed deviations. 

In the case of the oil shale * coke,” the figures 
for D, relating to a 700°C furnace temperature 
are in close agreement with each other, although 
they are much lower than for coal shale. A 
considerable increase appears at 800°C; this 
must be ascribed to the decomposition of the 
carbonates, present in relatively large quantities. 
Table IV shows that if the decomposition of the 
carbonates begins only at 700°, this phenomenon 
is nearly complete at 800° in the conditions of our 
experiments, i.e. with a constantly renewed 
atmosphere. Obviously the decomposition of the 
carbonates deeply affects the structure of the 
material considerably increasing the porosity. 


Table content in oil shale residue. 


Originally 
Burned at 700°C : 
after 30 minutes 7-32 
after 1 hour 7-29 
after 2 hours 6-93 
Burned at 800°C ; 
after 30 minutes 0-92 
after 1 hour traces 


Finally, it may be noticed that the values 
obtained for the diffusivity through the porous 
layer are considerably lower than the values for 
the usual gas diffusion film. This confirms that 
in the present conditions chemical resistance 
may be neglected in comparison with diffusional 


Combustion rate of coal shale 


207 


resistance, and that diffusion through the porous 
layer is effectively controlling. 

It is tempting to try to establish a relation 
between diffusivity and porosity ; some work has 
already been done in this direction by other 
authors, and the present author is engaged on 
this subject. The results of a tentative approach 
will be published in a subsequent paper. 
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NOMENCLATURE 
194 p T° 


a : 


b= [1 [1 + 2 
C = total carbon content, ©, weight 
C’ = carbon content from carbonates, °, weight of 


-10°3 


shale 
C” = COg content of carbonates, °, weight of shale 
Cy — combustible carbon content in original material, 


°, weight 
¢ = combustible carbon content at time @, °, weight 
c’ = combustible carbon content at time @, related to 
original material, °,, weight 


C, = specific combustion rate, g em?. sec 
D = diffusion constant for porous material, em? /seo 
Dy = diffusion constant for porous material by T°K, 


sec 
Dy = diffusion constant for porous material by 0°C, 
em? ‘see 
2 = energy of activation, cal mol 
k =a constant; = 1k 
L layer thickness of porous material, em 
Ny — molar rate of consumption of oxygen, mol ‘sec. em? 
P = Pog OXygen partial pressure in gas phase, atm 
oxygen partial pressure at carbon surface, atm 
py = film or layer pressure, atm 
P = total pressure, atm 


P, = weight of a lump of shale, g 
r — radius of sphere containing the unburned carbon, 
em 
ry = initial radius (or equivalent radius) of a lump of 
shale, em 
R = gas constant 
T = combustion temperature, K 
« = excess of air 
p = initial density of shale, g/em* 
pe — density at time 


= time, sec 
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Nature of wall heat transfer coefficient in packed beds 


Tuomas J. HANRATTRY 
University of Illinois, Urbana 


(Received August 1954) 


Summary—A mechanism of heat transfer from a tube wall to fluid flowing in a packed bed is 


proposed whereby : (1) a mass of fluid is thrown against the wall ; (2) the fluid at the wall assumes 
the wall temperature and transfer occurs inward ; (3) this mass of fluid moves away from the 
wall and is replaced by new fluid. On the basis of this mechanism the following equation was 
derived using order of magnitude considerations : 


k k 


This expression is compared with data on wall film coefficients for cylindrical packings. 
Good agreement is obtained. However, no such agreement is obtained for spherical packing. 
The explanation for this difference proposed by PLaurz [3] is subscribed to. 


Résumé— Pour le transfert de la chaleur entre la paroi d'un tube et un fluide en mouvement, 


dans un lit poreux, l’auteur fait les hypothéses suivantes : 
(1°) Une masse de fluide rencontre la paroi. 


(2°) Le fluide, au contact de la paroi, prend sa température et un transfert s’établit vers 


lintérieur. 


(3°) Cette masse de fluide quitte la paroi et est remplacée par une autre. 
Sur la base de ce mécanisme et en tenant compte de ces considérations par ordre d’importance, 


il obtient ’équation ci-dessus qu'il compare aux valeurs expérimentales des coefficients de * film 
le long de la paroi * pour des garnissages cylindriques. L’accord est bon, mais pour des garnissages 


sphériques, l'accord est moins satisfaisant. 
proposée par PaLutz. 


INTRODUCTION 


The interest in recent years in the development 
of an analytical method for the design of catalytic 
reactors has prompted many measurements of 
the rate of heat transfer from surrounding walls 
to a fluid flowing through a packed bed. These 
experiments have shown that for cases in which 
the D,/D, ratio is large enough such that packing 
non-uniformities do not exist the diffusion of 
heat within the bed can be described by a Fick’s 
Law type equation 


—> 

where K is an eddy thermal conductivity. How- 
ever, at the bounding wall a sharp break, resembl- 
ing a discontinuity, occurs in the temperature 


L’auteur accepte l’explication de cette différence 


profile. (See Fig. 1). Coperty and 
[1] have described this wall resistance in terms 
of a film coefficient : 


(2) 


Measurements of h for such a system with air 
as the fluid have been carried out by Fexrx [2] 
using Celite spheres and cylinders, and by PLautz 
[3] using glass spheres. The data for these two 
investigators for spherical packing are plotted 
in Fig. 2. Since these data represent a variation 
of tube diameter from 3 inches to 8 inches and a 
variation of particle diameter from } inch to 
} inch, the agreement is striking. 

The film coefficient varies with G7? and the 
proximity of the form of the equation describing 
these results to that obtained for heat transfer 


h(t, t;) 


Tuomas J. Hawnrarry : 


to a turbulent boundary layer suggests that the 
same type mechanism might be operative in the 
packed bed wall heat transfer. 


RADIAL DISTANCE, INCHES 


1. Temperature profile in a packed bed heat 
exchanger [1). 


Fig. 


Treatment of the data 


Nature of wall heat transfer coefficient in packed beds 


transfer occurs inward; (3) this mass of fluid 
moves away from the wall and is replaced by new 
fluid. Thus the heat transfer process is character- 
ized by a continuous replenishment of fluid at the 
wall caused by its deviation in flowing around the 
packing. By relating the rate of renewal of sur- 
face to the packing size and fluid velocity an 
expression for h was obtained which showed a 
dependence on 

On the basis of this, PLaurz suggested that 
the difference in the results obtained with spheres 
and cylinders arises from the fact that spheres 
possess only point contacts with the wall whereas 
cylinders offer a much greater obstruction to the 
flow. Therefore, in the case of spherical packing, 
channelling at the wall is a distinct possibility, 
while the flow adjacent to the wall with cylindrical 
packing might be described in the manner out- 
lined above. 

The possibility that this theory might clarify 
the divergence in wall heat transfer data suggested 
that its application to Friix’s data might be of 
interest. This has been done in this paper. 


Nature or THe Watt Resistance 


The break in the temperature profile cited as 
evidence for the existence of a wall resistance 


on cylinders yields two 
differences: (1) there is 
a decided dependence of 
AD upon the D, D, ratio; 


k 
(2) the film coefficient 
varies with G°*. Since the 
results with spherical and 
cylindrical packing were 
obtained in the range of 
Reynolds Numbers and 
D,, ratios some explan- 
ation is needed. 
Hanrarry [4] proposed 
a mechanism of wall heat 


T 


transfer whereby: (1) a 
mass of fluid is thrown 
against the wall; (2) the 
fluid at the wall assumes 
the wall temperature and 


Fig. 2. Wall coefficient data for packed beds of spherica, particles [2) [3). 
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might result from a variation in AK due to a 
variation in bed properties or from a slower 
rate of heat transfer from the wall to the system 
than in the bed itself. This paper discusses the 
latter cause. 

A radial variation of packing properties could 
be caused by a “ jamming” or “ bridging” of 
the particles. This effect could only be overcome 
if the tube were inserted through the packing 
in the same manner as a “cork borer.” The 
magnitude of this effect will vary with the D, D, 
ratio. At large values of D, D, the non uni- 
formities should be confined to a relatively small 
region in the vicinity of the wall. 

In fact, the data presented by Morr [5] on the 
variation of the fraction voids with D, D,, show 
that for spherical particles the fraction voids is 
independent of D; D,, at values greater than 10. 
However, for irregular particles Morr has shown 
that the variation in bed properties becomes 
unimportant only at much larger values of 
D, D,. 

For cases in which “ bridging” is negligible 
the “resistance to heat transfer will depend 
only upon the nature of the process by which 
heat is transferred to or from the wall. Two 
paths of travel are available, the solids or the 
fluid. Ferwix [2] has shown that use of packing 
materials of higher conductivity will increase 
markedly the wall coeflicient, h; so it appears 
that the solids might under certain circumstances 
be an important avenue of heat transfer. How- 
ever, for materials of low thermal conductivity 
the major portion of the heat interchange will 
occur through the fluid. 

This transfer process might be viewed in the 
same manner as in a turbulent boundary layer. 
Channelling would occur adjacent to the wall, 
and transfer in the layer results from turbulence 
generated at the tube and solids boundaries 
rather than by mixing of jets, as in the bed 
proper. 

However, another mode of transfer is suggestive 
if the fluid adjacent to wall does not flow in a 
continuous fashion but is “ sidetracked” by 
packing which might be in its path. The process 
might be viewed in the same manner as suggested 
by Danckwerts [6] in his analysis of film coetti- 


cients in mass transfer. Fluid of relatively 
uniform temperature is thrown against the wall 
and the fluid in this mass immediately adjacent 
to the boundary assumes the wall temperature. 
Transfer of heat occurs inward in the same manner 
as if a block of material were placed in contact 
with a hot wall. If the contact time of the fluid 
mass with the wall is very small then the heat 
will penetrate only a short distance inward. 
If the fluid in the immediate vicinity of the 
wall be unaffected by turbulence which might 
exist in the fluid mass, then transfer of heat 
will occur by molecular conduction. According 
to such a mechanism the wall resistance would 
depend upon the molecular properties of the 
fluid and the average contact time of these 
fluid masses with the surface of the tube. 

For a case in which * bridging ” effects are not 
existent this average contact time for a given 
type packing should be a unique function of the 
size of packing and the fluid velocity since it is 
the arrangement of the packing adjacent to the 
wall which would be responsible for the contin- 
uous movement of masses of fluid to and from 
the wall. An approximation of the average 
contact time and therefore of the film coefficient 
is possible by postulating this mechanism. 


MATHEMATICAL TREATMENT OF 

* Noncontinuvous” Liqurip 
The mathematical treatment of such a case 
developed by Danxckwerts [6] will now be 

presented, 

The rate of heat transfer into a fluid of infinite 
depth by molecular conduction would be described 

by the following solution of Fick’s Law : 


= (p Cite — Cyty) @) 


If an effective contact time of the fluid masses 
with the wall, @,, can be assumed, then the 
fraction of the surface with ages between @ and 
60+ dé@ is d0/0,. The rate of heat transfer 
at the column wall becomes 


de 
4 
(4) 


substituting for (@) 


ll 


2 1 k 
Ve Vc dx — 
If the wall film coeflicient be defined as 
a=h(t, 


h=2 


The average contact time may be approximated 
in the following manner: Assume the number 
of diversions of fluid flowing adjacent to the 
wall to be directly proportional to the tube 
length and inversely proportional to the particle 
diameter. 


then 


(5) 


L 


The constant C should be of the order of 
magnitude of one. The average contact time of 
the fluid in the system will be equal to L/U. 
Therefore the contact time per diversion or the 
average contact time of the fluid masses with the 
wall would be 


(6) 


LU D 
= 
Ct/D, CU 
Substituting into Equation (5) 
cU 
Vol J¢ (8) 


Rearranging in terms of dimensionless groups 


If C be assumed equal to one 
AD, D,Up  [Cyp 
k J k 


Wat Corrricrent Data 
The results of Piaurz [3] and of [2] for 
glass Celite spheres for air have been plotted in 
Fig. 2. The equation representing this graph is 

AD, = o12 por 

The results of Fe.rx on Celite cylinders have 
been plotted in Figs. 3, 4, 5, 6. These data are 


compared with Equation 10 by assuming U = U,/e 
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3 4 
Fig. 3. Wall coefficient data for 4” cylinders in 5” tube. 


and that = — 0.74. Equation 10 then becomes 


AD, = 0.95 D,Uoe 
k N 


A summary of this comparison is contained in 
Table 1. In all cases the data can be satis- 
factorily represented by a function of the type 
shown above. The numerical value of the 
estimated film coeflicient was within 50 per cent 
of the measured value. In fact the data for 
4 inch cylinders in a 5 inch tube and the data 
for } inch and j inch = } inch cylinders in a 
3 inch tube are in exact agreement. 


10% 2 3456 81 2 3 
Dp Us 
és 


Fig. 4. Wall coefficient data for j’ cylinders in 5°’ tube, 
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a 
102 2 4 6 86 103 2 4 6 
Op PU, 
Fig. 5. Wall coeflicient data for 3” = 4 cylinders in 5” 
tube. 


All of these results are for the heating of air 
flowing in an upward direction. Friix made 
some runs on the cooling of air and found that 
larger heat transfer rates were obtained than in 
the heating runs. This discrepancy is probably 
due to the existence of natural convection 
currents in the layers close to the wall for the 
cooling runs. Since no account is taken of this 
in the proposed mechanism, none of these data 
have been presented here for purpose of com- 
parison. 


100, 
— 
EQN. 10 
CYLINDERS 
102 234 68” 2 3 


Fig. 6. Wall coefficient data for 3” tube. 


Table I. Comparison of equation with experimental results. 


Estimated relation : 


AD 
= p~ oP 
k eu 


System Experimental equation 

AD Uop 

cylinders in 5” tube — 0-95 oP 
k 

DU op 

cylinders in 5° tube = 


eylinders in 5” tube 


cylinders in 3” tube 
= eylinders in 3” tube 


Discussion 

The agreement between experiments with cylinders 
and the wall film coefficients estimated by order 
of magnitude approximation using the “throwing” 
theory is much better than should be expected. 
In fact, the exact agreement in some cases may 
be fortuitous. It is entirely possible that some 
other mechanism might describe the results 
equally well. This possibility is not contended 
However, the fact that this theory gives a satis- 
factory interpretation of present experimental 
results justifies its consideration. 

Thus, as suggested by PLaurz [3], two types 
of boundary layers may exist adjacent to the 
tube wall, “ continuous” and discontinuous.” 
By making this assumption one can account for 
the marked difference in the wall coefficients 
measured for cylindrical and spherical packing 
produces discontinuous film by breaking up the 
fluid flow adjacent to the wall. Since spheres 
possess only point contacts with the wall this 
tendency is not so great. 

For spherical packing the Nusselt Group was 
found to be a unique function of Reynolds 
Number and this function appears to be inde- 
pendent of tube diameter or D,/D,, ratio in the 
range investigated. (D,/D,, = 12 to 40). 

However as shown by Fe.ix the tube diameter 
has a marked effect on the results with cylindrical 
packing in this same range of D,/D,, ratios. 
This is in agreement with the results of Morr [5] 
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who showed that a constant bed voidage is obtain- G = fluid mass velocity = py. 
ed at much larger Dy D, ratios for irregular h = wall resistance coefficient. 


ki This variation of ti Dine aaaen. k = fluid molecular thermal conductivity. 


ties of the bed might explain why the relation L = length of packed section. 
between Nusselt Group and Reynolds Number n — number of diversions. 
was found to vary with the packing size. da = heat transfer per unit area of tube wall. 
= temperature. 
ACKNOWLEDGMENT te = temperature of the wall. 
= fluid temperature at the wall as obtained from an 
The author wishes to express his gratitude to extrapolated temperature profile. 
R. H. Wine. for suggestions obtained from U = fluid velocity. 
him in discussions of the material contained in Uy = average fluid velocity based on empty tube 
this paper. cross-sectional area, 
N 6 time. 
average contact time of fluid masses with the wall. 
C, = heat capacity of the fluid. pw = fluid viscosity. 
C = constant. p = fluid density. 
Dy = tube diameter. (0) = rate of heat transfer per unit area from surface 
D, = particle diameter. of age 6. 
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Approximation methods for designing cooler-condensers 


R. C. Carrns 


School of Chemical Engineering, N.S.W. University of Technology, High Street, Kensington, Sydney, Australia. 


Summary~— The need for a general approximation method for the design of cooler-condensers 
is discussed and the recent papers on the condensation of vapour in the presence of a non-condensa- 
ble gas are reviewed. The various approximation methods that can be used for design are 
presented including the CoLpurn [5] method, in which enthalpy is assumed as the driving force. 

A new approximation method is given for obtaining the mean value of the heat flux, (U At), 
which is based on the assumption that a plot of the heat flux point values throughout the 
cooler-condenser, against the heat transferred per unit time, q, is a parabola. 

The areas calculated by the use of the various approximation methods are compared with 
the areas obtained for cooler-condensers designed by the rigorous CoLpuRN and HovucEn [7] 
method, for six examples that have been fully worked out in the literatare. 

It is concluded that the method involving the mean value of (U At) is considerably simpler 
to use than the CoLpurN and HouGen method and is also generally more accurate than the 
COLBURN approximation method for the cases studied. The error in the use of the CoLpuRN method 
varied from — 0-7°%, to + 95%. The error in the use of the mean heat flux method varied from 
— 07% to + 20°, except for one case which gave an error of — 42%. 

The methods involving various logarithmic means of the terminal values are shown to lead 
to serious errors. 


Résumé— L‘auteur insiste sur le besoin d'une méthode générale d’approximation pour le calcul 
des réfrigérants-condenseurs. Il passe en revue les publications récentes sur la condensation 
dune vapeur en présence d'inertes. Il présente les méthodes d’approximation valables pour ce 
calcul, y compris la méthode de CoLBURN dans laquelle l'enthalpie est prise comme force motrice. 
Il donne une nouvelle approximation pour obtenir la valeur moyenne du flux de chaleur (U At), 
approximation qui assimile 4 une parabole la courbe des valeurs ponctuelles du flux de chaleur 
d'un bout 4 l'autre du réfrigérant-condenseur en fonction de al chaleur transmise par unité de 
temps, q. 

Il compare les surfaces calculées en se servant des diverses méthodes d’approximation et les 
surfaces obtenues en appliquant la méthode rigoureuse de CoLBURN et HouGen, dans le cas de 
six exemples déja complétement étudiés dans la littérature. 

En conclusion, la méthode utilisant la valeur moyenne de (U At) est beaucoup plus simple 
que celle de CoLBuRN et HouGen et, généralement, plus précise que la méthode d’approximation 
de CoLtpurN dans les cas étudiés. L’erreur commise en utilisant la méthode de CoLBURN varie 
de — 0,7% a + 95%. Avec la méthode du flux de chaleur moyen, elle varie de — 0,7°% a + 2%, 
excepté pour un exemple qui donne une erreur de — 42%, 

L’auteur montre que les méthodes introduisant diverses moyennes logarithmiques des valeurs 
aux extrémités conduisent & de graves erreurs. 


INTRODUCTION 

The several papers appearing in the literature in 
recent years on the condensation of vapours 
from mixtures with non-condensable gases has 
reflected the growing interest in this subject. 

The method of Corpurn and Hovcen [7] 
has generally been accepted as the basis for the 
rigorous design of cooler-condensers to remove 
vapour from a gas-vapour mixture. 


This method has the disadvantage of being 
tedious in that it involves successive trial and 
error substitutions and because of this there is 
a need for an investigation of approximation 
methods. A method involving preferably only 
the terminal conditions which at least would 
always give a reasonable conservative design 
area is greatly desired. 

The nature of the problem makes the 
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attainment of an approximation method difficult 
since in a cooler-condenser both the temperature 
drop from cooling medium to the gas-vapour 
mixture and the overall coeflicient may vary 
widely throughout the condensation process. 
The stepwise procedure of CoLBuRN and 
allows for the variation of the temperature drop 
and overall coefficient by evaluating their product 
at several points in the cooler-condenser from 
equation (1) and graphically or otherwise integrat- 
ing equation (2), 


(UAt) = (t, —t.) + KM, A(p, — p.) = 
hg (t, — ty) 


dq 
| o (UAt) 

Sarru and Rowson [10] working with a water 
jacketed vertical condenser tube 10-23 feet 
long and 1-231 inch. LD., condensed benzene 
from a superheated benzene-air mixture for 
both upward and downward flow of the mixture, 
with the cooling water countercurrent in the 
jacket. The degree of superheat was from 18°F 
to 27°F. They applied the CoLsurn and Hovucex 
method of calculation to their experimental 
results assuming that the benzene-air mixture 
was at its saturation temperature. 

In the case of turbulent flow throughout 
the condenser and downward flow, the CoLBuRN 
and Hovcen method gave an area within 6°, 
of the actual area used. For other cases where 
turbulent flow was not maintained throughout 
condensation, the CoLpurn and Houcesx method 
gave conservative results when downward flow 
was used. An oversize area was also obtained 
for all but one case of upward flow. 

Cotpurn [5] has recently introduced an 
approximate procedure which is an application 
to indirect contact equipment of MerKe’s [9] 
total heat method for direct contact equipment. 
The method presupposes that enthalpy is the 
driving force for the transfer of heat. It is 
much simpler to use than the rigorous method 
since it involves only the terminal conditions. 

COLBURN writes :— 


— 
a= (U,/¢,) (Adan (3) 
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Where (U, is evaluated as the arithmetic 
mean of the two terminal values. The terminal 
values of (U, are found from 


c,/U, = ¢,/h, + (4) 
where m = di* dt 


= slope of the enthalpy saturation 
line. 


CoLBuURN obtained very good agreement between 
this approximate method and the CoLpurn and 
Hovcen method. 

In order to apply the and Hovcen 
method to the cooling of an unsaturated mixture, 
it is necessary firstly, to evaluate the operating 
line, as vapour is removed from the mixture. 
This in itself is a difficult procedure. Bras [1] 
has used the relationship, 


dp _ P, AP, ay (5) 
dt pam At, \Se_ 
previously proposed by CoLBuRN [5], to evaluate 
the operating line for the cooling of an unsaturated 
mixture of converter gases and water vapour in 
an ammonia oxidation unit. He evaluates the 
slope of the pressure-temperature operating line 
at the inlet of the condenser from equation [5] 
and follows the slope for a short distance. The 
slope is again evaluated at the new point and 
followed again. By continuing this procedure he 
finds the operating line on the pressure-tempera- 
ture curve. The stepwise method of CoLBURN and 
Hovucen is applied simultaneously for condensa- 
tion occuring along this derived operating line. 
Cairns [2] has shown that for the systems 
chloroform-air and steam-air which were studied, 
the Cotsurn and Hovcen method agrees to 
within 5°, of the area used experimentially. 
The gas-vapour mixtures were introduced very 
close to saturation and turbulent flow existed 
throughout the condensation process. The water 
jacketed condenser tube used was 6 feet long and 
1-081 in. I.D., with cooling water in the jacket 
countercurrent to the downward flow of the 
mixtures. CArRNs also compared the approxi- 
mate method of CoLBurRN to the systems and 
found that the areas were approximately 90°, 
greater than used experimentally. 
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In all of the above cases the sensible heat 
carried by the diffusing vapour has been neglected 
in the calculations. Co_tpurn and Drew [6] 
give a theoretical method for estimating the 
increase in sensible heat transfer obtained, when 
simultaneous heat and mass are transferred 
in the same direction, as in a cooler-condenser. 
These have not been the sensible 
heat transfer is only a small fraction of the 
total heat transfer. Further, Carrxs and Roper 
[3] have shown, in an experimental study of a 
wetted wall tower, that the theoretical method 
of CoLsurn and Drew cannot explain the data 
for the of simultaneous heat and 
transfer in opposite directions. 

In the absence of further experimental work 
on simultaneous heat and mass transfer in the 
same direction the comparisons made in_ this 
paper of the available approximation methods 
do not include the effect of the sensible heat 
of the vapour transferred on the heat transfer 
rates. 

There have been a sufficient number of examples 
of the application of the Co-nurn and HovuGen 
method, and various modifications of it, to 
enable comparisons to be made with this method 
and the several approximation methods that 
are available. The six examples chosen in this 
work cover a considerable range of physical 
properties of gas-vapour mixtures under a 
similarly wide range of terminal conditions. 


used since 


case mass 


APPROXIMATION METHODS 


The problem of designing cooler-condensers has 
resolved itself into two main Firstly 
is the method of CoLBuRN and HovuGen generally 
applicable and secondly is there an approximation 
method which can be used safely ? 

On the basis of present experimental results 
[2], [10], the and Hovcen method is 
accurate to within 10°, for fully developed 
turbulent flow at least. Hence this method 
ean be used as a standard of comparison with 
the approximation methods. 

The simplest and hence the most usual approxi- 
mation method is to evaluate the terminal overall 
coefficients and use the arithmetic average overall 
coefficient in conjunction with the logarithmic 


issues. 
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mean temperature difference. The accuracy 
of this method decreases as the difference between 
inlet and outlet overall coeflicients increases. 
Also the use of a logarithmic mean temperature 
difference theoretically requires that both the 
specific heat of the mixture and the overall 
coeflicient remain constant. However in cases 
where the coeflicient and specific heats vary only 
slightly this method is of value. 

The approximation method of CoLtpurn has 
been shown by CoLBurRN to give good agreement 
with the CoLpurN and Hovucen method in certain 
cases. It has also been shown to give as much 
as 90°, oversize in other cases. 

In some cases the use of the following method 
can be quite helpful. If the heat flux, (U At), 
can be shown or can be predicted to be linear 
or approximately linear with the heat transferred, 
q. then the logarithmic mean (U At) can be used 
to evaluate A from equation (6), 


= = 
(UAt)in 


(6) 

A disadvantage of this method is that, unless 
the conditions for which the method is applicable 
are approximately met, i.e. (UAt) linear in q, 
then either an undersize or oversize area will be 
obtained depending on the curvature. 

Another method, which suggests itself, is to 
evaluate the logarithmic mean of the products 
of the driving force at one end and the overall 
coetlicient at the other. This involves the assump- 
tion that the overall coefficient is linear in the 
temperature drop, as well as the temperature 
drop being linear in q. 

The foregoing approximation methods involve 
only the use of terminal conditions and 
thus simplify the design procedure considerably. 
However, none of the methods are generally 
applicable under all conditions since each has 
its own definite limitations, which are very often 
exceeded. 

The most promising results appear to be in a 
general method for obtaining a mean value of 
the heat flux, (UAt),, and using this to find 
the area. A mean value of (UAft) is desired 
since as previously stated both U and At vary 
throughout the cooler-condenser. 
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In searching for a general approximation 
method it has been noticed that when the (U At) 
point values are plotted against q the curve so 
obtained approaches very closely part of a 
parabola in shape, as in Fig. 1 for example, for 


p24 
= 
= 
= 
x 
= 
8 + 
| | 
4 
| | 
|| | 
° 2 4 8 
Fig. 1. (U At) vs. q for the Cotpurn and Hovucen [7] 


example. 


the Co.surn and Hovcen problem. This means 
that the relationship between (U At) and q can 
be represented approximately by a second order 
equation of the form, 


(U At) = + bq +e (7) 


On this basis it is possible to evaluate (U Af), 
to use in equation (8), 


qr 


A= (U Ad), (8) 


From equations (2) and (8), 


Subsitituting (7), 


(9) 


(10) 
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then from (8), or from equations (2) and (7), 
dq 
(aq? + bq + e) 
Integration of equation (11) gives ; 
Vb? \2aq +b + — 
when 6? > 4c, (12) 


A= (11) 


when b? < 4 ac, (13) 
—2 
and A = 7+ when b? = 4 ae. (14) 


When b+? > ac the extrapolated parabola of 
(U At) versus g would cut the axis at two places. 
When b? < 4ac the parabola does not cross the 
q axis. When b? = 4ac there is only one value 
of q when (U At) is made equal to zero. 

The solution of these integrals is simple, 
once values of a, b and ¢ are found. This may 
be done as follows :— 


Let (UAt) = (UAt),, at — dar 
(UAt) = (UAt) nig at = 
(UAt) = (UAt), at q= + 


Where (UAt),, = The heat flux at the 
gas-vapour outlet, 


(UAt)niq <= the heat flux at a point in 
the cooler-condenser where 
half the total heat has 
beer transferred and 

(UAt), = the heat flux at the gas- 


vapour inlet. 


Then from equation (7) 


When ¢ = 0 

= (UM) aia (15) 
When q = — 347 
(UAt), = 4 — + (18) 
When = + 347 
(Ut), = + + (UM)mig (17) 
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(17)-(16) ; 
(UA), —(UAt), = bay 
b = [(UAt), — (UAt),] Gr 
(17) + (16) ; 
(UAt), + (UAt),, = + 2 
a = 2[(UAt), + (UAt), — 2 (UAt)pia] (19) 


For these values of a, b and c¢, equation (11) 


A= 


(aq* + bq + ¢) 


~ 
2. 4 = In —— 
for b? > 4ac. 
a(_@+b ) 


“a7 


(21) 


Hence knowing the values of (U At) at each 
end of the cooler-condenser and a third value 
at a point in the cooler-condenser when half 
the total heat removed has been transferred, it is 
possible to calculate the values of the constants 
a, b and ¢ in equation (20). The appropriate 
integrated equation is then used, depending 
on whether b? = sac, which gives directly the 
area the conditions specified. 

The only limitation on this method is that 
(U At) when plotted against approaches a part 
of a parabola, 

Carey and Wiiiiamson [4] have recently 
published a graph, but without a mathematical 
development, for obtaining their mean driving 
force in direct contact air-water systems. This 
graph gives values of a factor, f with which to 
multiply the midpoint driving force in order 
to obtain the mean driving force, but it is limited 
to those values of the mean driving force which 


are greater than about half and less than 1-65 
times the mid value. The graph is also based 
on the assumption that the relationship between 
their driving force and the amount transferred 
is a parabola. 

In order to give a comparison of the approxima- 
tion methods mentioned here, they have been 
compared with the CoLtsurn and Hovucen 
method for those cases that have been fully 
worked out in the literature. 


APPLICATION OF APPROXIMATION METHODS 


Table I gives the conditions known for the 
examples that have been chosen from the litera- 
ture for this comparison. All but one case, 
the Bras example, are for the removal of vapour 
from saturated mixtures where the pressure- 
temperature operating line is immediately known. 

Table II gives the results of the application of 
the CoLBURN approximation method which is 
based on enthalpy as the driving force. The 
value of m for each case, except that of Bras, 
was evaluated at the temperature of the gas 
vapour mixture, ¢,, at inlet and outlet, by the 
method outlined previously [2]. The Smiru 
and Rosson example for the benzene-air mixture 
is shown, and in this example, it was necessary 
to correct the heat transfer coefficient at the 
gas-vapour outlet since the Scumipt number 
here was greater than the PranptL number. 

For the case of the unsaturated mixture of 
Bras, where water vapour is removed from 
ammonia converter gases, difficulties arise as to 
the value of m to use. In this example m at 
the top of the cooler-condenser cannot be 
evaluated as the slope of the saturation line 
at t,, as suggested by CoLpurn. It is therefore 
necessary to assign a position at which to evaluate 
m, which was taken as the point (t,, i,*). 

The correct slope to use will actually be greater 
still than that at (t,, i,,*) but if the area obtained 
in using the slope at (¢,, i,*) is greater than 
the area obtained by the CoLBuRN and HouGEen 
method, any increase in m will over design the 
cooler-condenser further. As is shown later an 
oversize area is obtained by the Co_BuRN 
approximation method even when m is evaluated 
at (t,, i,,*) for the top of the cooler-condenser. 
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becomes 
dq 
—| 
for b? < 4ac. (22) 
and A = a. when b? = 4ac (23) 
|| 
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Table 1. 
and and Bras CAIRNS CAIRNS KERN 
Hovucen Rosson example example example example 
example example No. 1 No. 2 
Saturated Saturated Unsaturated Saturated Saturated Saturated 
System nitrogen-steam air-benzene converter air-chloroform air-steam carbon-diovide- 
vapour gas-steam vapour steam 
Gas-vapour In 203-0 169-7 278-0 120-2 168-3 267-0 
te °F Out 104-0 90-3 85-0 94-0 129-9 120-0 
Cooling 
water In 7740 85-6 67-0 58-5 66-4 80-0 
temp. — — 
Out 140-0 60-3 70-3 115-0 
Gas rate Ib. 
moles hr. 100-0 0-0576 246-68 0-20 0-869 35 
Cooling water 
rate Ib. hr. 154,000 2760 172,000 2010 2010 137,000 
Total pressure atm. 1 1 1 1 1 3-05 
Exchanger type Shellandtube, Single tube Shellandtube, Single tube Single tube — Shell and tube, 
1-1. Gas on vertical. 1-1. Gas on vertical. vertical. 1-4. Gas on 
shell side. tube side. shell side. 
Vertical 
Gas-vapour flow Countercurrent Downwards Downwards Downwards Downwards 
and and and and 


countercurrent countercurrent countercurrent countercurrent 


The value of m at the gas-vapour outlet in 
this case although unsaturated was taken as 
the slopes of the enthalpy curve at (, = 85°F, 
as before, which is safe. 

The recalculated values of Cairns for chloroform 
and air, and steam and air are also shown in 
Table II together with the example of Kern [8] 
on the removal of steam from carbon dioxide. 
In this case, since the condensation is carried 
out under pressure, allowance is made for the 
absolute pressure in evaluating m from the 
equation previously presented [2]. 

In Table II] the results are shown of the 
application of the method based on the assump- 
tion that a plot of (U At) versus q is a parabola, 
in which the area is obtained directly by integra- 
tion. The values of (UAt) at each end of the 


cooler-condensers were taken from the references, 
as was the mid-point for the saturated cases. 
This midpoint value can be obtained of course 
without knowing all the (UAft) point values 
for these cases. 

Figs. 1 to 6 show the relationship between the 
calculated point values of (U'At) and q for each 
of the examples. Fig. 3 plots the point values 
of Bras using the (U Af) values as calculated 
by Bras. 

The deviation from the straight line in the 
chloroform-air example of Carrns is due to the 
exclusion of the condensate resistance at the 
point of the gas-vapour inlet. The dotted 
line indicates approximately the relationship 
between (U At) and qg at this end of the cooler- 
condenser. 
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Bras CAIRNS CAIRNS KERN 
No. 2 


Smiru and 
Rosson 


| Top Bottom | Top Bottom Top Bottom Top Bottom Top Bottom 


850 1202 940 1683 1299 2670 | 1200 


Vapour pressure | 
Py atm. 0-883 0-178 0-151 0-038) 0-664 O-151 2-68 O-115 
| | | 


Enthalpy, i, above 


32°F 4787) 136-8 193-3 40-1070 485 147 8556 36-3 

(12640 «1182 105 101 | 26606 235 

Ai 4661186 1638S HS 232 458 123 3523 21-2 


66-1 12:3 17-7 413 494 0-72 


h, corrected for Se. 


(U,/¢,) arith. avg. 


Area by and 
Hovucen method ‘ 


Error 


Fig. 6, for Kexn’s example, illustrates clearly the curve for this example does not approach 
the existence of a point of inflection. As is seen that of a parabola. 
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Table 2. 
| | 
| 
Gas-vapour temp. 
Cooling water temp. | | | | | 
ty °F | 880 | 85-6 | 670 | 58-5 | 703 664 115-0 80-0 
L. 
54 
| | — | | 
fe | | | 
h, 481 00694 5-37 302 B13 574 189 | 57 
he 150 113 113 179 98 159 159 926 926 
0-48 | 108 27-7 27-2 1-60 308 61 13-1 1-62 | 25-4 
B.t.u. /hour 7139 «1,077,400 3781 7795 4,802,000 
0-78 37-5 oe 18-5 
840 63-9 515 255 685 
Calculated area A | | 
sq. ft. 10-9 613 314 | 3-18 519 
—= 
| 
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Btu /(hrlisq ft) 


VAt —e 


6 
Btu/hr 


Fig. 2. (U At) vs.q for the Surru and Rowson [10] example. 


For the midpoint value in the Bras example, 
where the operating line is initially unknown, 
it is necessary to make the simplifying assumption 
that the operating line between inlet and outlet 
points is a straight line on the humidity chart. 


7 


Btu ft) 


2 3 a 
Btu/he 


(UAt) vs. for the Carrns [2] chloroform-air 
example. 
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BtuAheXsq ft) 


UAt 


4 6 Le] 12 
Btu/he 


Fig. 3. (U At) vs. q for the Bras [1] example. 


It is also necessary to assume that the outlet 
is the same composition and temperature as 
was assumed by Bras. For the purpose of 
comparing the methods on the same basis this 
second assumption is justified, although where 
the outlet concentration is not known for the 
cooling of an unsaturated mixture, it would be 
necessary to assume the mixture left saturated. 

The temperature and humidity at the point 
in the cooler-condenser where half the total 


10 


Btu/(hriisq ft) 


UAt x 10° 


° a 
Btu/hr 


Fig. 5. (U At) vs. ¢ for the Cams [2] steam-air example. 
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Table Ill 


and | Ssuru and 
| Rosson 


Bras 


CAIRNS 
No. 1 


Gas-vapour 
outlet (U At); 


Use equation 


| 

| 
Gas-vapour | | 
inlet (UAt), ‘18,540 9704 4295 5030 8120 97,950 
(U 23,500 4032 2650 “2000 "5380 "54,000 
ar) B.t.u./hr. 9,706,000 "1,077,400 9795 "4,802,000 
a — 5696 x 109-10 6525 x 9-35 x 1-161 x 10°* — 9875 x 10°° 8-15 x 10°! 
x 10 1.066 2,026 x 10-2 
3.035 x ie 1-839 1-206 x 0-551 105 x 
4ac | x105| 1052 | 1207 | — 1-76 x 10% 


Calculated area 


heat has been transferred were calculated by an 
enthalpy balance, neglecting the sensible heat in 


UAt x 


° 2 5 $ 
Btu/hr 


Fig. 6. (UAt) vs. ¢ for the Kern [8] example. 
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the condensate, and from the assumed linear 
relationship between H and ¢,. The value of 
(UAt),ig Was then found by trial and error 
substitution in equation (1). 

The value of (U Af)... obtained in this way was 
2650 as compared to the value obtained by Bras, 
by rigorous calculation, of 2720. The tempera- 
ture of the mixture was found to be 184°F as 
compared to 175°F by Bras. 

The application of this method of evaluating 
the area can best be shown by the following 
calculation of the CoLsurnx and HouGeN example: 


Values previously calculated by CoLpurn 
and Hovucen ; 
(UAt), = 1,629 B.t.u./(hr.) (sq. ft.) 
(UAt), = 18,540 B.t.u. (hr.) (sq. ft.) 
(UAt),.iq = 23,500 B.t.u. (hr.) sq. ft.) 
= 9,706,000 B.t.u. /(hr.) 


from equations (19), (18) and (15) ; 


| 
| | CarRns Kern 
mis 
IL. 
3 
(2!) (21) (22) (21) (21) 
A sq. ft. | 690 5-92 | 572 | 170 | 220 | 
Error | — 07% | — 05% + 2.0% | + | + | — 42% 
| 
3 a 
2s 
|| 
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a = 2(18,540 + 1629 — 2 x 28,500) / From Fig. 8 reference (4), 
(9-706 = 10°)? f = 0-595 
= — 5-696 x (UAt),, = 0-595 238,500 
b = (18,540 — 1629) (9-706 » 10°) = 13,980 
- 17-42 x 10 A = 9,706,000 13,980 
= 2-35 = 694 sq. ft. 


Table IV gives the results of the method 
using the graph of Carey and WILLiamson, 
which enables solution of the appropriate integral 
= — 5.353 x 10°% to be read off graphically. 

Table V compares the various other methods 
involving logarithmic means of the terminal 


and 4ac = — 4 x 5.696 « 10° » 
2.35 « 10° 


*. use equation (21) 


— tac = 75-21 values only. 
b — 4/b® — 4ac = — 57-79 x 10% In Tables II to V the percent error is expressed 
b + 4/68 — sac = + 92-63 x 10+ as the difference between the areas obtained by 


the Conpurn and Hovcen method and the 
= approximation methods, divided by the area 
75-21 x 10% obtained from the and Houcen method, 
— 5-696 x 10° g — 57-79 « 10°*\] *er all multiplied by 100, 
5-696 x 10° x + 92-63 x 10 
= 690 sq. ft. ComPaRIsSON OF APPROXIMATION METHODS 


A 


Table IL shows that the CoLBeurNn approximation 
Using the chart presented by Caney and gives an oversize area varying from 9-3°%, to 


WILLIAMSON, 95°,. The area obtained by CoL_surNn using 
UAt his approximation method for the CoLsurn 

( 
(Uat), 144 and HovuGen example was 690 sq. ft., an error 
of — 0.7%, the Cotpurn and HovuGen method 
(UM) nia — 1.97 giving 695 sq. ft. On the basis of the examples 
(U compared here, the approximation 


Table 


| 
COLBURN and and Bras CAIRNS CAIRNS KERN 
HovuGen Rosson Nov 1 No. 2 


(UAt)mia At), 


(UAt)mia/(U 1-27 0-42 0-62 0-52 0-66 0-55 


Factor, f, Fig. 8 [4] 0-595 Not on graph O-7T15 0-855 O-88 Not on graph 


-condensers 
— 8-085 1076 
VOL. 
3 
a 1954 
fr 144 175 5-74 2-60 2-30 82-9 
(UM), «13,980 1,895 2023 | 4738 | = 
A sq. ft. 694 _ 369 | 1-70 1-65 — 
Error O1% — + + | + | 
| 
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Table V 
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Bras 


and Rosson 


COLBURN 
and HouGen 


Top Bottom 


Area by 


1290 (+ 86% )* 


4-29 (— 28%) 571 (+ 1-8°%,) 
mean of 
(Ug At,) and 
(Uy 


1733 (+ 149%) | 15-1 (4+ 154%) | 536 (—4-5%) 


Area by 


(UMD 1396 (4 101%) 446 (— 25%) 627 (4 12%) 


CAIRNS 


84-0 


CAIRNS 
No. 2 


59-9 


35-5 


1-45 (— 14%) 


98-0 


63-5 152-0 


1-63 (0-0%) 173 (— 54%) 


1-66 (— 1-2%) 


1-52 (— 


1-78 (+ 92%) 


483 ( +27%) 


1-68 (4+ 31%) 247 (— 85%) 


* Percentage in brackets is the percent error based on the area obtained by the CotpurRN and HovuGeEn calculation. 


method appears to yield on oversize area rather 
than undersize and is thus conservative from a 
design viewpoint. However a cooler-condenser 
with 95°, more area than is required is an 
expensively designed —cooler-condenser. As 
COLBURN points out the method is only intended 
to give a rough approximation for the true area. 

In Table III the results of the method based 
on the parabolic relationship between (U At) 
and q is seen to give areas from — 0.7%, + 2-0°%, 
of the Cotsurn and Hovucenx method in all 
cases except the Kern example, which gives a 
— 42°, error. Apart from this exception, the 
agreement is remarkable considering the reduction 
in the volume of calculations afforded by this 
approximation method as compared to the 
CoLpurN and Hovcen method. 

The small error of only + 2°, for the Bras 
example supports the assumptions made in order 
to apply the approximation method. If the actual 
figure for the mid-point value of (U At) of 2720 
is used the exact area is obtained, which would 
be expected by inspection of Fig. 3. 

The reason the method gives a 42°, undersize 
area for the example of Kern is due to the non- 
parabolic relationship between (UAt) and 4, 
as has been previously mentioned. This in 


turn is due to the plot of the heat flux due to 
latent heat transfer against g giving a point of 
inflection. In all other examples this quantity 
shows the normal parabolic relationship. It is 
difficult at this stage to predict when the non- 
parabolic relationship would occur, since only 
one out of the six examples shows this tendency. 

Table IV gives the areas as obtained by the 
graph of Carey and WILLIAMSON for those cases 
where the graph can be used. The agreement 
with Table III is good. The Carry and 
WILLIAMSON graph is recommended for the cases 
where the ratios of the mid-point (U At) values 
to the terminal (UAt) value can be plotted. 
This procedure eliminates the necessity of solving 
equations (21), (22) or (23) and reduces further 
the amount of calculations. Where the Cargy 
and WILLIAMSON graph cannot be used the area 
can be obtained from equations (21), (22) or (23). 

The results of the methods involving the various 
logarithmic means of the terminal values are 
shown in Table V. The correlation of areas by 
these methods is shown to be generally poor. 
The use of these methods must be treated with 
considerable caution as either an unreasonably 
undersize or oversize area can be obtained. 

It is seen from Tables Il and III that the 
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| | | 
No. 1 
U 118-8 4-91 20-6 30-8 28-2 82-9 36-9 O44 
at 630 270 
| 
| 
| 
34 
|| 
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method using (U' At), is generally better than 
the CoLBURN approximation method, except 
for the example of Kerx. It must be emphasised 
however that three point values of (U' At) are 
required involving three trial and error subsitu- 
tions, whereas the Co_nurN method does not 
require any trial and error substitutions and 
operates entirely on the terminal conditions. 

The use of only three point values of (U Af) 
in the mean heat flux method however means 
a considerable reduction in the amount of calcula- 
tion compared to the usual six to ten point 
values. 

It must be mentioned that the evidence 
offered here is not conclusive but would have to 
be augmented by further worked examples before 
it could be definitely shown that the mean heat 
flux method is more suitable in all cases than 
the CoLBURN approximation method. the 
basis of the present work it appears to be so in 
five out of the six cases. 

An attempt was made to modify the CoLpurn 
approximation method by obtaining a mean 
value of the product of the enthalpy heat transfer 
coeflicient and the enthalpy driving force, 
(Ue, Ad cane in the same manner as (U Af)... 
This was based on the assumption that this 
product showed a parabolic relationship with q. 
In the case of the Cornurn and Houcen example 
the area obtained was 3,240 sq. ft. compared to 
695 sq. ft. actual, and consequently this method 
was rejected, 


CONCLUSION 

An approximation method has been presented 
for the design of cooler-condensers. This method 
is based on the assumption that a plot of the 
point values of (U At) against q, is a parabola, 
The method has been shown to give good agree- 
ment with the Cotpurn and Hovcen method 
in all the cases studied except one. The non- 
agreement in this case is discussed. 

In applying the method the graph of Carry 
and WILLiamMson is recommended. When this 
graph cannot be used the appropriate integrated 
equation can be used which gives the area 
directly. 

This method is less tedious to use than the 


Cotsurn and Hovcen method and it is more 
accurate in five out of the six cases studied 
than the Co_surn method of approximation. 
More examples need to be solved however before 
its limitations could be fully determined. 
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NOTATION 


a = constant defined by equation (19). 

A = cooler-condenser surface area, sq. ft. 

b = constant defined by equation (18). 

¢ = constant defined by equation (15). 

= humid heat, B.t.u. /(Ib.) 

dq = inerement of the total heat transferred, 
B.t.u. hour. 

h = heat transfer coefficient, B.t.u./(hr.) (sq. 
ft.) (°F), 
hy — the combined conductances other than 

the gas film, 

h, = the gas film. 

il = humidity of the gas-vapour mixture, 

Ib. dry gas). 

enthalpy of the gas vapour mixture, 

B.t.u. /Ib. 

i.* < the saturation enthalpy at ¢,, 

i, = the enthalpy at ¢,, 

i,* = the saturation enthalpy at ¢,,, 
(Ai), = the overall logarithmic mean 
enthalpy driving force. 

K = mass transfer coefficient, 
(sq. ft.) (atm.). 
m = slope of the enthalpy saturation line. 
== di* /dt. 
M,, = molecular weight of the vapour. 
p = partial pressure, atm., 
Ppm ~ the logarithmic mean of p, and 
(P-p,). 
— the vapour pressure at 
Pg — the non-condensable gas in the main 
stream, 
Pp, = the vapour pressure in the main 
stream, 
Ap, = the partial pressure driving force 
across the gas film. 
P = total pressure, atm. 
Pr = Prandtl No. 
q = heat transferred per unit time, B.t.u./hour, 


I! 


226 


1954 


R. C. Camns: Approximation methods for designing cooler-condensers 


qr = total heat transferred. 
Sc = Schmidt No. 
t = temperature, °F, 
t. = condensate surface, 
ty = gas-vapoyr mixture, 
to = water, 
At = overall temperature driving force. 
At, = temperature driving force across the 
gas film. 
U = overall heat transfer coeflicient, B.t.u. 
/(hr.) (sq. ft.) (CF). 


= arithmetic average of U, and U4. 


(U,/e,) = overall heat transfer coefficient per unit 
enthalpy difference, B.t.u./(hr.) (sq. ft.) 
(B.t.u. /Ib.). 
(UAt) = heat flux B.t.u./(hr.) (sq.ft.) 
(UAt),, = the true mean (U At). 
(UAt),,, = the logarithmic mean of the (U At) 
terminal values, 
(UAt)mia = (UC At) at the point when half the 
total heat has been transferred. 
A = latent heat of vaporisation, B.t.u./Ib. 
Subscripts : 1 = gas-vapour outlet, 
2 = gas-vapour inlet. 
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Book review 


J. M. Covtson and J. F. Ricwarpson : Chemical En- 
gineering, Volume I, 1954, Pergamon Press Ltd., London 
and McGraw-Hill, New York. VIII + 370 pages. Price 
38s. 6d. ($7.50). 


Tue main purpose of this book is to serve as a text for 
chemical engineering students and, in the reviewer's 
opinion, it will be extremely useful in this purpose. Apart 
from that, it is gratifying that an European counterpart 
of similar American textbooks, such as G. G Brown's 
“ Unit operations,” is now available. 

The entire book will consist of two volumes, in which 
the fundamentals of the physical operations in chemical 
engineering are treated. In view of this, the title * Chemical 
Engineering * suggests a somewhat wider scope than is 
actually covered in the book. Volume I deals with fluid 
flow and heat and mass transfer. In Volume II these 
foundations will be applied to the design of equipment for 
physical separation. 

More precisely, Volume I is divided in four sections : 
A - Flow of fluids, B — Heat transfer, C - Momentum, heat 
and mass transfer, and D — Humidification and water 
cooling. These sections are preceded by a useful chapter 
on “ Units and Dimensions” in which the C.G.S., the 


F.P.S., and the British engineering systems are com- 
pared. 


Personally, I regret the absence of a discussion 
of the rationalized Giorgi-system (based on kg-mass, metre, 
second, degree Centigrade and Ampére), since this system 
is frequently used on the continent. 

Section A contains a chapter on the energy relationships 
in flow systems and one on friction in pipes and channels. 
These topics are treated in the conventional way, but the 
typical hydraulic applications of flow in open ducts are 
relatively more stressed than in other books of this kind. 
However, the reader will not find a treatment of the 
important group of phenomena which include the force 
exerted by a fluid flow on particles, the settling of spheres 
and otherwise shaped particles, flow through packed beds, 
ete. The authors have, apparently, preferred to introduce 
these subjects in the second volume. The next chapter 
treats the various means of measuring flow rates in 
industrial plants, and is followed by a chapter on pumps 
for chemical works. As the authors have pointed out, the 
latter is somewhat outside the scope of the book but it 
has been inserted because pumps play such an important 
role in all flow systems. Unlike the others, this chapter 


is primarily descriptive in nature, but elementary con- 
structional aspects have generally been omitted. 

Section B is entirely devoted to heat transfer. Due 
attention is given to all the basic mechanisms, with 
special emphasis put on the application to the design of 
heat exchangers. In the section on unsteady heat con- 
duction, however, it would perhaps have been better if 
the operational calculus for solving Fourier’s equation 
had been left out in favour of a number of important 
results, e.g., the well-known graphs for heating or cooling 
of slabs, cylinders and spheres where the dimensionless 
Fourier number is used as a parameter. 

Section C is attractive because of its original approach. 
After a very clear chapter on diffusion and the basic 
concepts of mass transfer a comparison is made between 
momentum, heat and mass transfer. It is to be regretted 
that the latter chapter contains only 14 pages, whereas 
a number of applications, especially in the field of mass 
transfer, might have been treated in order to give the 
students more opportunity to grasp the importance of 
the material in this chapter. The next chapter, which 
treats the boundary laver, is more of theoretical interest ; 
the boundary laver theories do not yet provide a solid 
basis for chemical e. .ineering calculations. Among other 
things, the “ universal velocity profile * of turbulent flow 
along a solid boundary still seems to be open to criticism. 

The final section of Volume I deals with simultaneous 
heat and mass transfer for the system water-air. The 
essential difference between the wet-bulb temperature and 
the adiabatic saturation temperature has been clearly 
pointed out. The section ends with a discussion of the 
methods of making calculations for cooling towers, which 
is very good, but which could be improved by the inclusion 
of values of the transfer coeflicients for a few typical 
cases. 

Throughout the text the reader is referred to carefully 
chosen literature references for further study. Many 
illustrative examples, but no problems are included. The 
paragraphs containing more advanced theoretical treat- 
ment are printed in small type. 

Summarizing, a thorough text-book, carefully written 
and well printed, with rather much theory and mathe- 
matical treatment, at the same time a good introduction 
to several important design problems. 


Ii. Kramers (Delft) 


